Prineil  Services  Technical  Information  Agency ' 

Because  of  our  limited  supply,  you  are  requested  to  return  this  copy  WHEN  IT  HAS  SERVED 
YOUR  PURPOSE  so  that  it  may  be  made  available  to  other  requesters.  Your  cooperation 
will  be  appreciated. 


NOTICE:  WHEN  GOVERNMENT  OR  OTHER  DRAWINGS,  SPECIFICATIONS  OR  OTHER  DATA 
ARE  USED  FOR  ANY  PURPOSE  OTHER  THAN  IN  CONNECTION  WITH  A DEFINITELY  RELATED 
GOVERNMENT  PROCUREMENT  OPERATION,  THE  U.  S.  GOVERNMENT  THEREBY  INCURS 
NO  RESPONSIBILITY,  NOR  ANY  OBLIGATION  WHATSOEVER;  AND  THE  FACT  THAT  THE 
GOVERNMENT  MAY  HAVE  FORMULATED,  FURNISHED,  OR  IN  ANY  WAY  SUPPLIED  THE 
SATD  DRAWINGS,  SPECIFICATIONS,  OR  OTHER  DATA  IS  NOT  TO  BE  REGARDED  BY 
IMPLICATION  OR  OTHERWISE  AS  IN  ANY  MANNER  LICENSING  THE  HOLDER  OR  ANY  OTHER 
PERSON  OR  CORPORATION,  OR  CONVEYING  ANY  RIGHTS  OR  PERMISSION  TO  MANUFACTURE, 
USE  OR  SELL  ANY  PATENTED  INVENTION  THAT  MAY  IN  AITY  WAY  BE  RELATED  THERETO. 


Reproduced  by 

DOCUMENT  SERVICE  CENTER 

KNOTT  BUILDING.  DAYTON,  2.  OHIO 


THE  DYNAMIC  RESPONSE  OF 
FLOATING  BRIDGES  TO  TRANSIENT  LOADS 


By 


J . P.  Romualdi  - Project  Engineer. 

E.  O'AppoIonia  - Associate  Professor 

F.  T.  Mavis  - Head 


Department  of  Civil  Engineering 
Carnegie  Institute  of  Technology 


For 


Office  of  Naval  Research 


Contract  No.  N onr  - 760(03) 


August,  1954 


1 


THE  DYNAMIC  RESPONSE  OF  FLOATING  BRIDGES  TO  TRANSIENT  LOADS 

By 

J.  P.  Roniualdi  - Project  Engineer 

E.  D'Appolonia  - Associate  Professor 

F.  T.  Mavis  - Head 

Department  of  Civil  Engineering 
Carnegie  Institute  of  Technology 


August  195^ 


• « 

! 

\ 

I 

H This  work  was  sponsored  by  the  Army,  Navy  and  Air 

^ I Force  through  the  Joint  Services  Advisory  Committee  for 

Research  Groups  in  Applied  Mathematics  and  Statistics  by 
^ Contract  No-  Nonr-760(03 ) - 

. .. 


i 


TABLE  OF  CONTENTS 


Page 

ABSTRACT vi 

ACKNOWLEDGEMENTS vii 

INTRODUCTION 1 

PREVIOUS  WORK 7 

NUMERICAL  SOLUTIONS 11 

General 11 

Single  Ponton  Bridge 13 

Bridges  of  More  Than  One  Floating  Support 20 

Damping 26 

Accuracy  of  Numerical  Procedure 27 

Numerical  Example 31 

ELECTRONIC  ANALOG  SOLUTIONS 32 

Results  of  Analog  Solutions 34 

Accuracy  of  Analog  Solutions 39 

IMPULSE-DISPLACEMENT  TESTS 40 

SUMMARY 45 

CONCLUSIONS 49 

APPENDIX  A 51 

APPENDIX  B 55 

APPENDIX  C 58 

BIBLIOGRAPHY 66 

VITA  67 


LIST  OF  FIGUHES 


Figure  No. 

2. 

3. 

4. 

5. 

6. 


Sub.lect 

Schematic  Diagram  of  Ponton  Bridge 
Simple  Beam  Analog 
Single  Mass  Oecillator 
Single  Ponton  Bridge 

Buoyant  Force -Displacement  Relationship  for 
an  M-4  Ponton 

Reaction  on  Ponton  due  to  Moving  Load  P 


7.  Displacement-Time  Relationship  for  Single 
Ponton  Bridge 

8,  Virtual  Mass-Displacement  Relationship  for 
an  M-4  Ponton 


9.  Displacement-Time  Relationship  for  Single 

Ponton  Bridge  Considering  Virtual  Mass 

10,  Relationship  Between  Impact  Factor  and  Time 
of  Crossing  for  Single  Ponton  Bridge 

11,  Relationship  Between  Impact  Factor  and  Time 
of  Crossing  for  Single  Ponton  Bridge  Consi- 
dering Mass  of  the  Moving  Load 

llA.  Relationship  Between  Impact  Factor  and  Time 

of  Crossing  for  Single  Ponton  Bridge  Consi- 
dering Mass  of  the  load 

12,  Two  Ponton  Bridge 

13,  Displacement-Time  Relationship  for  Two  Ponton 

Bridge 

14,  Displacement-Time  Relationship  for  Three 

Ponton  Bridge 

15,  Displacement-Time  Relationship  for  Three 

Ponton  Bridge  Considering  Virtual  Mass 


ii 


LIST  OF  FIGURES  (continued) 


16. 

17. 

18. 

19. 

20, 

21. 

22. 

23. 

2A. 

25. 
25A. 

26. 

27. 

28. 

29. 
29A. 

30. 


Displacement-Time  Relationship  for  Three 
Ponton  Bridge  Considering  Virtual  Mass  and 
Damping 

Numerical  Integration 

Displacement-Time  Relationship  for  Single 
Mass  Oscillator 

Influence  of  Moving  Load  on  Ponton 

Displacement-Time  Relationship  for  Single 
Ponton  Bridge  Considering  Mass  of  the  Load 

Displacement-Time  Relationship  for  Single 
Ponton  Bridge  Considering  Mass  of  the  Load 

Displacement-Time  Relationship  for  Single 
Ponton  Bridge  Considering  Virtual  Mass  and 
Mass  of  the  Load 

Influence  of  Moving  Load  on  Pontons  1 and  2 

Displacement-Time  Relationship  for  Two  Ponton 
Bridge 

Displacement-Time  Relationship  for  Two  Ponton 
Bridge  Considering  Virtual  Mass 

Relationship  Between  Impact  Factor  ard  Time 
of  Crossing  for  Two  Ponton  Bridge 

Displacement-Time  Relationship  for  Two  Ponton 
Bridge  Considering  Virtual  Mass  and  Damping, 
C“  100  Ib-sec/in. 

Displacement-Time  Relationship  for  Two  Ponton 
Bridge  Considering  Virtual  Mass  and  Damping, 
C*=  loco  Ib-sec/in. 

Influence  of  Moving  Load  on  Pontons  1,  2, 
and  3 

Displacement-Time  Relationship  for  Three 
Ponton  Bridge,  ElA>^  “ 85-8  Ib/in. 

Relationship  Between  Impact  Factor  and 
Time  Crossing  for  Two  Ponton  Bridge 

Displacement-Time  Rf '' ationship  for  Three 
Ponton  Bridge,  El/L-^  *=  42.9  Ib/in. 


iii 


LIST 

31. 

32. 

33. 

34. 

35. 

36. 

37. 

38. 

39. 

40. 

41. 

42. 

43. 

44. 

45. 

46. 

47. 

48 . 

49. 

50. 

51. 


OF  FIGURES  (continued) 

Displacement-Time  Relationsldp  for  Three 
Ponton  Bridge  Considering  Virtual  Hass  and 
Damping,  C“  100  Ib-sec/in. 

Disp].acement-Time  Relationship  for  Three 
Ponton  Bridge  Considering  Virtual  Mass  and 
Damping,  C“  1000  Ib-sec/in. 

Displacement-Time  Relationship  for  Five 
Ronton  Bridge 

Comparison  of  Numerical,  Exact,  and  Analog 
Solutions  for  Two  Ponton  Bridge 

Virtual  Mass  Apparatus 

Schematic  Diagram  of  Virtual  Mass  Apparatus 
Test  Pontons 

Schematic  Diagram  of  Teat  Pontons 

Buoyant  Force-Displacement  Relationship  for 
Teat  Pontons 

Calibration  Curve  of  Test  Beam 

Virtual  Mass-Displacement  Relationship  for 
Test  Pontons 

Model  of  M-4  Military  Ponton 

Impulse-Displacement  Apparatus 

Schematic  Diagram  of  Impulse-Deflection 
Apparatus 

Results  of  Impulse-Displacement  Test  On 
Rectangular  Ponton 

Results  of  Impulse-Displacement  Test  on 
Rectangular  Ponton 

Results  of  Impulse-Displacoment  Test  on 
Triangular  Ponton 

Results  of  Impulse-Displacement  Test  on 
Triang-ular  Ponton 

Results  of  Impulse-Displacement  Teat  on 
Triangular  Ponton 

Results  of  Impulse-Displacement  Test  on 
Semi-Circular  Ponton 

Results  of  Impulse -Displacement  Test  on 
Seiai-Circul  ar  Ponton 


iv 


LIST  OF  FIGURES  (continued) 

52.  Results  of  Impulse-Displacement  Teat  on 
Semi-Circular  Ponton 

53.  Relationship  Between  Impact  Factor  and  Number 
of  Pontons 

54 • First  and  Second  Modes  of  Vibration  in  One, 

Two,  Three  and  Five  Ponton  Bridges 

55 • Relationship  Between  Impact  Factor  and  Time 

of  Crossing  for  Single  Ponton  Bridge  Consi- 
dering Mass  of  the  Lead 

56.  Displacement-Time  Relationship  for  Single 
Ponton  Bridge  Considering  Mass  of  the  Load 

57.  Relationship  Between  Impact  Factor  and  Ratio 
of  K to  EI/L3 

58.  Relationship  Between  Impact  Factor  and  Time 
of  Crossing  for  Two  Ponton  Bridge  - C*0,  100 
and  1000  Ib-sec/in. 

59.  Relationship  Between  Impact  Factor  and  Time 
of  Crossing  for  Three  Ponton  Bridge  - C“0, 

100  and  1000  Ib-sec/in, 

60.  Relationship  Between  Impact  Factor  and  Damping 
for  Two  and  Three  Ponton  Bridges 

61.  Determination  of  Damping  Coefficient 

62.  Basic  Computing  Components 

63.  Block  Representation  of  Relationship  Between 

and  Y2 

64.  Block  Representation  of  Relationship  Between 
Yi,  Y2  and  Y3 

65.  Solution  of  the  Differential  Equation 

+ BPyi  + Dy^  _ f(t)  ■=  0 

66.  Solution  of  the  Differential  Equation 

P2yi  t-  Byip2y^  + - f(t)  - 0 

67.  Programming  of  Two  Ponton  Bridge 


V 


ACKNOWIEDGMEflTS 


This  research  project  described  herein  was  sponsored  by  the  Army, 
Navy  and  Air  Force  through  the  Joint  Services  Advisory  Committee  for 
Research  Groups  in  Applied  Mathematics  and  Statistics  by  Contract 
No.  Nonr -760(03 ) . The  authors  acknowledge  their  indebtedness  to  the 
members  of  the  Civil  Engineering  Department  for  their  helpful  sug- 
gestions, and  to  Robert  Hagaman  and  Dayton  Cook  for  their  aid  in 
carrying  out  the  computations.  They  are  also  grateful  to  Dr.  H.  M. 
Trent  and  Dr.  W.  A.  McCool,  both  of  the  Naval  Research  Laboratories, 
Washington,  D.  C.,  for  their  kindness  in  making  available  the  analog 
computor  at  the  Naval  Research  Laboratories. 


Vi 


IMROOLTmON 


A floating  bridge  is  essentially  a continuous  bridge  on  n 
intenT'.ediate  supports  which  act  as  elastic  springs.  The  resistance 
of  a support,  brought  about  by  the  bu<^ant  force  of  the  water,  may 
be  linear  or  non-linear  with  respect  to  displacement  depending  upon 
the  geometry  of  the  support.  A study  of  the  c^ynamic  behavior  of 
such  structures  is  coinp-licated  by  the  fact  that  the  mass  of  a 
floating  support,  or  ponton,  is  not  constant.  This  is  due  to  the 
virtual  mass  of  the  j onton  which  must  be  considered  in  the  case  of 
bodies  accelerating  in  fluids.  The  virtual  mass  of  the  ^onton  is 
the  )iass  of  the  [.-onton  in  air  ,vlus  the  added  ,nuass  which  is  usually 
expressed  in  terns  of  the  i.ass  of  displaced  water.  The  added  mass 
arises  from  the  vertical  acceleration  of  the  ponton  through  the  water. 
However,  as  the  ponton  is  disi laced  the  quantity  of  displaced  water 
changes  and  hence,  the  added  mass  changes,  '..'e  are  then  confronted 
with  the  problem  of  a mass  which  varies  with  displacement.  In  the 
literature  Uiere  is  little  treatment  of  tliis  ty;.e  of  non-linear  pro- 
blem. By  foriiial  macheniatical  methods  solutions  of  engineering  pro- 
blems of  this  type  are  intractable,  f-xist  non-linear  problems  treat 
the  case  of  mass  varying  v;Lth  time,  as  in  the  | roblems  associated  with 
rocket  flight. 

The  m.ost  important  application  of  a dynamic  analysis  of  floating 
bridges,  at  present,  is  in  the  design  of  militaiy  i onuon  bridges.  A 
brief  survey  of  existing  design  standards  indicates  that  these  bridges 
are  designed  for  static  loads  only  whereas,  as  will  be  subsequently 
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shovin,  the  djmajiiic  dis  >lacai>ents  for  a ^±vtm  moviri^  ioi-d  can  be  Lvdce 
as  {treat  as  the  corresponci'i;^  static  displaceitiGiitl  This  is  im;X)rtant 
not  only  vdth  respect  to  tlie  maximum  y-ermissible  displacsmsi.t  before 
swamping  but  also  vdth  rtispect  to  the  stresses  set  up  in  the  ponton 
shell. 

This  investigation,  however,  is  concerned  vdth  displacements  and 
results  are  presented  in  terms  of  iin[>act  factor®.  The  impact  factor 
is  defined  as  the  ratio  of  the  maxim'Jm  djTiaiidc  displacenent  of  the 
bridge  due  to  a moving  load  to  the  maximum  static  displacement  that 
cam  be  caused  by  the  same  load  if  it  v/cre  apidicd  statically  to  the 
bridge.  For  bridges  vdth  an  odd  number  of  pontons  the  maximum  static 
displacement  occurs  viien  the  load  is  at  t- e center  ponton.  The  maximum 
d.Tiamic  displacement,  hovrevor,  usually  occurs  after  the  moving  load 
has  passed  the  center  ponton.  If  t)>e  L;ajdmum  stresses  dec  to  a moving 
load  are  desired,  they  can  be  calcuiatcd  to  a good  ai'proxLmation 
from  the  c^amic  defiected  load  line  by  determining  tlie  static  forces 
necessary  to  produce  the  same  displacements.  The  stresses  c:_n  then 
be  calculated  from  the  bending  moments  produced  by  Oiese  loads. 

A floating  ponton  bridge  is  a series  of  discrete  floating  units 
normally  vd.th  some  degree  of  articul.ation  between  the  declv  structure 
and  the  pontons.  Articulation  is  not  considered  in  this  work.  It  is 
assumed  that  the  bridge  is  floating  at  its  dead  v^eight  displacement 
which  is  5-2  in.  for  an  K-A  ponton  bridge.  /J.1  displacements  pre- 
sented in  this  vjork  are  measured  from  this  datum  and  the  structure 
is  considered  continuous.  This  means  f ’at  tin  resistance  to  dis..'lace- 
ments  due  to  moving  loads  cones  from  three  sources;  (1)  the  flexural 


rigidity  of  the  continuous  dec'.,  structure,  (2)  Uie  t force  of 

the  \>Qter  end  (3)  ibe  inertia  of  tlie  structure.  The  buoycjit  reslstunce 
of  a ponton  to  displacement  is  15CX)  lbs/ in.  'I’he  coal'ined  resistance 
to  displacement  of  a single  ponton  bridge  of  30  ft  span,  due  to  the 
buoyant  force  of  the  water  and  the  flexural  rigidity  of  ti;-;  bridge,  is 
34420  Ibs/in.  These  figures  give  the  relative  o.  1 :r  of  i.a^nitude 
of  the  resistance  of  the  deci,  structure  to  resistatce  of  the  water. 

The  capacity  of  a single  ponton  before  swamping  is  about  6C000  lbs 
which  is  about  eight  times  the  \Jeight  of  a ponton  plus  one-h.alf  the 
adjoining  spans.  The  mass  of  the  loads  sxijiported  by  a ponton  bridge, 
therefore,  is  largo  in  comparison  to  the  mass  of  the  bridge.  This 
study  has  been  extended  to  take  into  account  the  mass  of  the  load  for 
loads  up  to  78000  lbs. 

The  natural  period  of  vibration  for  i.ionton  bridges  ranges  from 
0.150  sec  for  a single  ponton  bridge  to  0.64  tec  for  a five  ponton 
bridge.  These  values  neglect  virtual  :ass.  The  actual  periods  are 
less  due  to  the  effect  of  virtual  mass  and  are  different  at  different 
depths  of  submergence.  Hovrever,  to  serve  as  a basis  of  comparison, 
the  tine  for  a moving  load  to  cross  a given  structure  is  given  in 
terms  of  the  natural  fundiimental  period  of  that  sti'ucture  neglecting 
virtual  mass.  In  general  the  ...aximum  dynariiic  effects  occur  Uien  the 
load  crosses  the  briige  in  about  one  or  tira  tires  the  fundajnental 
period  of  the  bridge.  For  a three  ponton  bridge,  for  examiile,  \d.th 
a natural  fundamental  [,’eid.od  of  0.37  sec  and  a span  length  of  bO  ft 
the  speed  of  the  moving  load  vd.ll  be  55  mph  if  the  load  crosses  in 
tvdee  the  fundamental  period.  For  a five  ponton  bridge  (socUi 
length  = 90  ft)  the  corresponding  speed  is  49  mph.  The  r:jige  of 
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times  of  crossing  studied  was  selected  to  include  this  nonnal  range 
of  operation. 

In  general,  a ponton  bridge  consists  of  many  lontons. 

But  the  analysis  of  bridges  of  more  than  four  or  five  pontons  is 
extremely  laborious.  The  number  of  simultaneous  differential  equations 
describing  the  beha’.-ior  of  such  structures  is  equal  to  the  number  of 
pontons.  This  work,  however,  is  concerned  only  with  bracketing  the 
effects  of  the  different  variablee that  affect  the  behavior  of  such 
floating  structures  and  tc  establish  trends  and  show  relative  variations . 
The  study  of  bridges  of  from  one  to  five  pontons  is  considered  suffi- 
cient for  this  purpose. 


A continuous  floating  bridge  on  n intermediate  supports,  sub- 
jected to  a load  F moving  across  the  span  vdth  velocity  v,  is  shown 
in  Fig. 
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FIG.  1 SCHKFATIC  DIAGRAM  OF  PONTON  BRIDGE 
The  general  equations  describing  the  displacement  of  each  sup- 
port, 1,  2,  3...n,  due  to  the  moving  load  I,  are 

a^y^  = ?^{t) 

^ ^21^1  " ^"22"K)72-^ a^y^  = F^(t)  (l) 


or 


V ^ (t) 

n nl  1 n2^  2 nn  n n 


n 

(jih  " i 


(2) 


i=l 
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where 


6^j,  the  Kronecker  Delta,  is  0 when  i/j 

and  1 when  i=j 

iTij^  is  the  equivalent  concentrated  mass  at  support  i, 
(ibs-sec^/in. ) 

is  the  reaction  at  support  i due  to  a unit  displace- 
ment at  support  j,  (Ibs/in.) 

K is  the  slope  of  the  buoyant  force-displacement  relation- 
ship of  the  floating  support,  (ibs/in.) 
is  the  displacement  of  support  i,  (in.) 

F^(t)  is  the  external  force  at  support  i due  to  the  load  P 
at  any  time  t,  (lbs) 

The  general  cases  considered  in  this  work  are: 

(a)  and  K are  constants.  The  load  P is  a constant 
force  (its  mass  is  neglected).  is  a linear  func- 
tion of  the  position  of  the  load  P on  the  span;  and 
since  F moves  across  the  s;an  v-dth  constant  velocity 
Fj^  is  linearly  time  dependant.  The  equations  then 
appear  as  in  (2). 

(b)  ra^'  is  the  virtual  mass  and  is  given  by 

^ y^ 

The  equations  of  motion  then  become 

n_ 


(c)  The  mass  of  the  moving  load  P is  considered.  Hence, 
is  a function  of  time  and  the  inertia  force  of  the 
moving  load,  F,  The  equations  of  motion  can  then  be 
represented  as 


n 


' -F^(t,yp|=  0 


(4) 


where  yp  is  the  ve  bical  acceleration  of  the  bridge 
at  the  position  of  the  load  P. 

(d)  Viscous  damping  is  considered.  Letting  C be  the  co- 
efficient of  viscous  damping  the  equations  are 


A =1 


(e)  The  most  general  case  includes  virtual  mass,  damping, 
and  the  mass  of  the  moving  load.  The  equations, 
assuming  a linear  spring  constant  K,  are 


n 


i=l 


In  this  thesis  a numerical  procedure  developed  by  J.  M.  Bonetti 
is  applied  to  the  solution  of  these  problems.  The  method  is  based  on 

the  concept  of  impulse -momentum  and  the  solution  is  carried  out  by 

successive  integration  over  short  time  intervals.  The  method  is 

extended  to  the  non-linear  class  of  problems  of  cases  (b),  (c),  and  (e). 
(1)  Numbers  in  raised  parentheses  refer  to  Bibliography 


Problems  in  each  of  the  above  five  categories  were  solved  by  this 
numerical  method  as  illustrations.  The  bulk  of  the  solutions,  however, 
were  obtained  on  the  electronic  analog  computer  which  is  installed 
at  the  Naval  Research  Laboratories  in  Washington,  D.  C.  The  method 
of  solution  by  the  analog  computer  is  discussed  and  in  a few  cases 
the  analog  and  manerical  solutions  are  compared.  Exact  solutions 
of  the  differential  equations  of  motion  by  formal  mathematical 
methods  are  obtained  to  check  the  analog  and  numerical  solutions. 

Experimental  tests  were  conducted  to  verify  the  importance  of 
the  influence  of  virtual  mass  on  the  transient  oscillations  of 
floating  str'jctures.  From  known  impulses  calculated  solutions  were 
obtained  neglecting  virtual  mass  and  damping  and  considering  virtual 
mass  and  damping.  These  results  are  compared  with  experimental 
results.  An  experimental  procedure  for  determining  virtual  mass  is 
also  presented. 


PREVIOUS  \i(ORK 

The  beha\dor  of  a floating  bridge  subjected  to  moving  loads  falls 
into  the  broader  classification  of  moving  loads  on  simple  span  bridges. 
Indeed,  a floating  bridge  can  be  thought  of  as  a simple  span  bridge 
whose  flexural  rigidity  is  augmented  by  the  addition  of  intermediate 
springs. 

Early  investigations  of  the  effect  of  moving  loads  on  bridges 
were  conducted  by  Sir  G.  G.  Stokes  who  was  concerned  with  the  deflec- 
tion of  a simply  supported  girder  subjected  to  a concentrated  mass 
mo-ving  \vith  constant  speed  Applications  of  his  treatment  are 
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limited  because  the  mass  of  the  load  is  assumed  largo  compared  to 
tho  mass  of  the  girder  which  is  neglected  in  his  method  of  solution 
The  case  of  a constant  force  moving  across  a simply  supported 
girder  was  studied  by  Sir  C.  E.  Inglis^^^  and  by  S.  Tomoshenko^^^ . 

In  this  case  the  differential  equation  of  motion  for  a simple  span, 
subjected  to  a moving  load  P whose  mass  is  neglected,  is 

_oL!*--Z.  7 = 1 P(x,t)  (7) 

^ t2  El 

where  w is  the  weight  per  unit  length  of  girder 
El  is  the  flexural  stiffness  of  the  girder 

= w/EIg 

y is  the  vertical  displacement  at  distance  x along  the  span. 

In  general,  the  moving  force  P(x,t)  can  be  expressed  as  an 
harmonic  series  whence,  equation  (?)  becomes 
h 2 ^ 

^ y + k^  Jl Z_  = 2P \ sin  iff  x sin  i ff  vt  (8) 

EH  L L 

i=l 

where  L is  the  length  of  the  span 

V if  the  velocity  of  the  moving  load  P 
The  general  case  in  which  the  mass  of  the  load  and  the  mass  of 
the  beam  are  considered  was  first  investigated  by  H.  H.  Jeffcott^^^. 
The  method,  however,  based  on  successive  approximations,  was  cumber- 
some and  did  not  converge  in  mar»y  Cases.  A more  complete  analysis, 
based  on  the  method  of  Inglis,  was  afforded  by  A.  Hillerborg^^\  He 
also  expressed  the  forcing  function  as  an  harmonic  series  and  obtained 
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the  differential  etjuations  of  motion  for  a simple  span  in  the  form 

3 r d^y(vt,t)  y sinJJUfl  sin  iwx  (9) 

^ EH  ^ Ji^  L L 

where  ql  is  the  mass  of  the  moving  load. 

This  is  essentially  the  same  as  the  equation  obtained  by  Inglis 
except  that  the  reaction  of  the  load  against  the  beam  is  changed  by 
the  amount 

„ _ d^y(vt,t) 

dx^ 

g 

This  is  due  to  the  effect  of  the  beam  "falling  away"  from  beneath 
the  moving  load. 

The  solution  of  problems  by  these  methods  is  laborious  and  is 
often  a needless  refinement  over  quicker  numerical  methods  which 
yield  approximate  solutions  to  sufficient  accuracy  for  engineering 
purposes.  One  of  these  methods  consists  of  representing  the  actual 
structure  by  an  analog  as  shown  in  Fig.  2. 


m 


m 


\ 



L 

r 

FIG.  ? SIMILE  PEAM  AGALOG 


The  mass  of  the  beam  is  assumed  concentrated  at  n intermediate 
supports.  The  remaining  portions  of  the  spa.i  are  considered  weight- 
less but  otherwise  retain  their  elastic  properties.  Considering  the 
external  force,  reaction  and  inertia  force  at  each  concentrated 
mass,  n equations  are  obtained  for  the  n intermediate  points  in  the  form 
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^ *12^2 


V2  "22^2 


a 


2n^n 


= Fi  (t) 

= (t) 


(10) 


“r/n  " Wl"  ‘^n2y2  ^ Wn  “ ^2 

The  greater  t'le  member  of  intermediate  points  n,  at  which  the 
mass  of  the  beam  is  assumed  concentrated,  the  closer  will  be  the 
approximation  to  the  actual  span. 

A convenient  stepwise  method  of  solution  for  equations  (lO)  was 
developed  by  J.  M.  Bonetti^^^  The  method  coneists,  essentially, 
of  assuming  the  impulse  at  each  intermediate  point  to  be  linear  ovei 
a short  time  interval  At  and  equating  this  impulse  to  the  change  in 
momentum  of  the  corresponding  concentrated  mass  at  each  intermediate 
point.  The  displacement  y^  of  any  intermediate  mass  at  time 
t+At  is  found  in  terms  of  its  displacement  and  velocity  at  time  t. 

The  extension  of  tJiis  method  to  floating  bridges  is  easily  per- 
formed. Equations  (lO)  are  the  same  as  equations  (l)  for  a floating 
bridge  with  the  exception  of  the  K values  of  equations  (l).  Moreover, 
the  analog  described  above,  which  is  only  an  approximation  of  a simple 
span  bridge,  is  actually  a much  closer  approximation  of  a floating 
bridge  since  part  of  the  mass  of  the  bridge  is  actually  concentrated 
at  the  intermediate  supports  due  to  the  presence  of  the  pontons.  The 
major  difference  lies  in  the  fact  that  the  concentrated  .masses  of 
the  floating  bridge  are  not  constant  but  vary  with  displacement. 

The  nvunerical  method  mentioned  above  has  been  extended  to  take  into 
account  these  varying  massos. 
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NUMERICAL  SOLUTIONS 


General 

The  numerical  procedure  used  is  based  upon  the  principle  of 
impulse  - momentum.  The  integration  of  the  resultant  force  on  a 
body  over  a time  interval  is  equal  to  the  net  change  in  the  momentxjm 
of  the  body.  In  this  case  the  integration  is  performed  numerically  - 
that  is,  the  integration  is  carried  out  successively  over  short  time 
intervals . 

Consider  the  mass  m suspended  from  a spring  vdth  spring  con- 
stant K shown  in  Fig.  3* 


FIG.  3 SINGIi:  M.\£5  CSCILD\TOR 

If  a force  P,  which  varies  ^nth  tine,  is  applied  in  the  positive 
(downward)  y direction  the  equation  o'f  motion  for  the  mass  is 

P (t)  - Ky  - m ^ (11) 

The  unbalanced  force  at  any  instant  is  the  difference  between  the 
applied  force  P and  the  tension  in  the  spring  Ky,  or 

UP  = P (t)  - Ky  (12) 

Introducing  this  expression  in  (11)  and  learranging  terms  we  have 

IF  dt  •=  mdv  (13) 

The  Integral  of  the  unbalanced  force  from  time  t^  to  time  t2  is 
equal  to  the  total  change  in  momentum  over  that  time  interval. 
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,to 


IF  dt  = m (V2-  =•  mAV 


(14) 


If  the  time  interval,  t2  - is  very  small  it  can  be  assumed 
that  the  unbalanced  force  varies  linearly  vdth  time.  With  this 
assumption,  letting  At  = t2  - t]^,  the  impulse  is 

'^2  ^ / ^1 
IF  dt  - 2 -1- 


(15) 


The  integral  of  the  velocity  of  the  mass  from  t^^  to  t2  is  equal 
to  the  change  in  displacement,  ^2  ~ yi>  durLng  the  time  interval. 

st.-> 


V d t = y.  - yn 


(16) 


Assuming  that  the  velocity  varies  linearly  »-dth  time  over  the 
interval  At  equation  (16)  becomes 
\ t.^ 


V d t 


= V2  * Vi 


At  = 


^2  - 


and  this  allows  the  change  in  velocity  AV  to  be  expressed  iri  the  form 


AV  = V2  - V^_ 


2 

It 


ivz  - - ^1^^) 


(17) 


Substituting  equations  (15)  and  (l?)  into  equation  (14)  the  impulse 
and  change  in  momentum  over  a short  time  interval  At  are  related  by 


IF  2^  IFi 


At  = (y^  _ y^  ^ V^At) 

At 


(IS) 


Using  expression  (12)  and  rearranging  terms  the  equation  for  the 
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displi'jcement  at  time  t2  is 


(JiS— +K)y2  “ (4m  - K)yi  V^At  + Pi  + ?2  (19) 

At^  At^  At2  ^ 


Thus  the  displacement  at  the  end  of  any  time  interval  is  expressed 
in  terms  of  the  displacement  and  velocity  at  the  beginning  of  the 
time  interval  and  the  applied  load. 

Single  Ponton  Bridge 

Fig.  4 represents  a single  ponton  bridge. 


L/2 


L/2 


FIG. 4 SINGLE  PONTON  BRIDGE 


The  spring  constant  K due  to  the  buoyant  force  of  the  displaced 
water  is  obtained  from  Fig.  5 which  shows  the  displaced  weight  of 
water  as  a function  of  displacement  for  a full  M-4  Ponton.  Fig.  5 is 

appended  at  the  end  of  this  thesis.  The  curve  is  not  linear  but  a 

value  for  K of  15CC  Ibs/in.  is  a good  approxijTiation . The  force  K'  to 
produce  unit  displacement  of  the  ponton  is  resisted  by  the  buoyant  fore 
and  the  flexural  rigidity  of  the  span.  For  a 30-foot  span  with  El  = 

320  X 10®  lb-in. ^ value  for  K'  is  34420  Ibs/in.  The  mass  of  the 

ponton  plus  one-half  the  mass  of  the  span  is  20.2  Ib-sec^/in.  Choosing 
the  time  interval  At  as  one-twentieth  of  the  fundamental  period  of  the 
bridge,  that  is 

At  = T„/20  = iE  / m_ 

0 10  V K' 

The  coefficient  of  4m/At^  becomes  400K'^^I? 
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Substituting  these  numerical  values  into  (19)  the  displacement  at 
time  t2>  neglecting  virtual  mass,  is 

* l)y.  “ ( ^ -l)y.  UOO  v,At  * ^1*^2 

7T2  ^ 7T2  Jf2  1 irr— 


(20) 


Tha  reaction  of  the  load  P on  the  ponton,  as  the  load  moves  across 
the  span,  is  assumed  to  vary  linearly.  The  beam  force  which  transmits 
the  effect  of  the  load  to  the  ponton  is  the  same  as  that  of  a span  com- 
pletely articulated  at  the  ponton.  In  other  words  continuity  is 
neglected  when  computing  the  force  at  a ponton  due  to  a moving  load. 

It  must  be  emphasized,  howe';er,  that  in  computing  the  resistance  to 
displacement  by  this  force  continuity  is  not  neglected.  Letting  H 
be  the  time  for  a moving  load  to  cross  the  span  the  reaction  of  the 
load  P on  the  ponton,  as  a function  of  time,  is  represented  by  Fig.  6 


The  results  for  a load  P = 23,400  lbs  (three  times  the  weight 
of  the  ponton  and  one-half  the  deck  structure  concentrated  at  the  center) 
for  different  speeds  of  crossing  are  shown  in  Fig,  7.  The  ordinate 
is  the  ratio  of  the  dynamic  displacement  to  the  static  displacement  that 
would  occur  if  the  load  P were  applied  statically  at  the  center  of  the 
span.  The  abscissa  is  time  in  multiples  of  the  funda:aental  period,  Tq. 
Example  calculations  for  H = T^  are  given  in  Appendix  A. 

If  the  virtual  mass  is  considered  equation  (19)  becomes 

+K')y2  = (..■^.  -K')yi  " JK-  K'ViAt  " Pi  " P2 

At2  At2  ^ 


(21) 


Keeping  Tq  as  the  time  scale  and  choosing  At  as  one-twentieth  the 
fundamental  period  Tq  equation  (21)  becomes 

(40.568  -HL  +I)y2  = (40.568  _E1  -l)yi  + 40.568  _El  V^At  + (22) 

m m m K' 

Fig.  8 shows  the  virtual  mass-displacement  relationship  for  an 
M-4  Ponton.  The  data  were  obtained  experimentally  from  a one-twelfth 
scale  model.  Since  added  mass  is  a function  of  shape,  mass  of  dis- 
placed fluid  and  direction  of  motion,  the  added  mass  of  the  full  scale 
ponton,  at  a given  di.splacement,  is  obtained  from  the  added  mass  of  the 
scale  model  at  the  corresponding  displacement.  The  shape  and  direction 
of  motion  is  the  same  in  each  case.  The  added  mass  of  th<j|inodel,  there- 
fore, is  multiplied  by  the  ratio  of  the  mass  of  displaced  water  of  the 
full  scale  ponton  at  the  given  displacement  to  the  mass  of  displaced 
water  of  the  model  at  the  corresponding  displacement  to  obtain  the  added 
mass  ol  the  full  scale  ponton. 

The  observed  added  mass  of  the  ponton  is  a linear  function  of 
displacement.  A more  detailed  description  of  the  metliod  of  obtaining 
these  values  if  given  in  a later  section. 

From  Fig.  8 the  virtual  mass  for  the  ponton  is 
m'  = 7A  + 2.23y 

when  the  initial  dead  v/eight  displacement  of  the  ponton  is  5.2  in. 
Choosing  ra  = 20.2* Ibs-sec^/in. , the  coefficient 

40.568  — = 148.613  - 4.478  y 
m 


Equation  (22)  becomes 
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(149.613  + 4.478y,)72  - (147.613  4.478y,)y, 

(23) 

♦ (148.613  + 4.478/1 ) ViAt  + 

K' 

It  is  assumed  that  the  virtual  mass  during  any  interval  of  time,  At, 
is  a function  only  of  the  displacement  at  the  beginning  of  the  time 
interval. 

Results  of  solutions  of  equation  (23)  are  given  in  Fig,  9 for 
H “ Tq,  for  a moving  load  P of  39>000  and  234>000  lbs  (five  and 
thirty  times  the  concentrated  weight  at  the  ponton,  respectively). 

An  example  solution  is  shown  in  Appendix:  A.  Fig.  10  summarizes 
solutions  of  equation  (23)  for  a moving  load  P of  three  times  the 
concentrated  weight  at  the  center  of  the  span  and  for  different  times 
of  crossing.  The  ordinate  is  the  impact  factor,  and  the  abscissa  is 
the  ratio  of  the  time  of  crossing  to  the  natural  period  of  the  structure 
ignoring  virtual  mass. 

If  the  mass  of  the  moving  load  is  considered  the  vertical  dis- 
placement of  the  bridge  at  the  point  of  application  of  the  moving  load 
must  be  known  in  terms  of  the  displacement  at  the  center  of  the  bridge. 
Assume  a linear  rel-j._on  between  the  deflection  y^^  at  any  position  x 
along  the  bridge  and  the  displacement  y of  the  ponton. 

whore  2x/L  when  0 < x < L/2 
and  o<  = 2(L-x)/L  when  L/2  <,  x <.  L 

If  the  load  is  assumed  to  follow  the  deflection  curve  of  the 

bridge,  the  acceleration  of  the  load,  when  0<x<L/2,  is 

• • • • . 

yx  ''K  y + V (24) 


-i  ■. 
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and  when  L/2  < x <.  L 

(24a) 

^ L 

where  v is  the  horizontal  velocity  of  the  moving  load  as  it  crosses  the 
bridge  and  V is  the  vertical  velocity  of  the  ponton.  The  last  term  in 
equations  (24)  and  (24a),  when  multiplied  by  the  mass  of  the  moving 
load,  is  the  Coriolis  force  vinich  is  due  to  the  precession  of  the 
velocity  vector  v as  the  span  displaces.  In  general  this  force  will 
not  have  much  effect  at  low  velocities  of  the  moving  load. 

The  force  of  the  moving  load  P applied  to  the  bridge  at  any  posi- 
tion X is  given  by 

assuming,  as  above,  that  the  force  F at  the  ponton  due  to  the  force 
varies  linearly  with  x 

F = F^o<  = (p  _ mLyx)«< 

The  unbalanced  force  on  the  ponton  is  the  difference  between  this 
downward  force  of  the  moving  load  and  the  upward  force  due  to  the 
flexural  rigidity  of  the  span  and  the  buoyant  force  of  the  water.  These 
latter  forces  depend  on  displacement. 

IF  = F - K'y  = (P  - mL/y)©^-  K'y  (25) 


If  the  Coriolis  force  in  equations  (24)  and  (24a)  is  neglected, 
equation  (2^  becomes 

IF  - Po<-  111^04^  y - K'y 

Substituting  this  expression  into  equation  (I4) 


f^2  2 

j (P«<-  K'y)dt  (m+mj^<»4 
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and  integrating  vre  have 

po<^  - K'y^  ^ * “H*^)(y2  - yi  - ViAt) 

SiiTiplifying,  the  expression  for  the  displacement  at  time  t^  is 


(Jjm_ 

. + Ai“L 

* K')y2  - 

( ta 

* - K')yi 

At^ 

At^ 

At^ 

At^ 

+ (4m  + 

Af2 

At 

4mL 

At^ 

o<  2)V^At  + P (o<^  +o<2) 

Letting  m 

At  - •»)  Tq  - iti'/f:’ 

2 2 2 2 2 2 2 
(1  73  n-  )y^  - (1  o<  - >]  TT  )y^  +(1  +y3d.)V3_At 

. 4?^  ^ K cat) 

Solutions  of  equation  (26)  are  shown  in  Fig.  11  fory3  “ 1»  3,  5,  and 
10  and  for  different  speeds  of  crossing.  An  example  solution  of  equa- 
tion (26)  is  given  in  Appendix  A. 

If  the  Coriolis  force  in  equations  (24)  and  (24a)  is  not  neglected 
equation  (25)  becomes 

IF  ■=  Po(  - m,  y*  - 4iJV  «><  V VA  - K'y 

for  0 <,  X < L/2  and 

IF  * P«»<  - y t 4nij^  <k  V V/L  - K'y 


for 


L/2  < x < L 
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juc;itituting  these  expressions  into  equation  ill,) 

r^2  2 

I (Po<,  - K'y)dt  = (m  )^V  ♦ , c>< ( V dt 

J ti 


for 

0 < X < L/2  and 

r 

2 

(P«4  - K'y)dt  •=  (ra  + 

h 

L/2<.x<L 

1 

for 

Integrating  and  rearranging  terms  wo  have 

Vdt 


2 _ 2 


+ (1  P (°<L  **^2)  ^27) 

for  0 < X < L/2  and 

? 2/3>)  ,.2  2 2 2/377,  2 2. 

^ )y2  * ^ )7l 

+ (1  )Vj_At  ■*•  rr  p ■*®^)  (27a) 

K ' 


whore  H =/<  T 

0 

Results  of  equations  (27)  and  (27a)  fory3  5 JUid  for  different  times 
of  crossing  are  shown  by  cui*ve  A in  Fig.  11a.  These  results  are  com- 
pared to  solutions  of  equations  (26)  which  neglects  the  Coriolis  force 
foryS  ■ 5 and  different  times  of  corssing  (curve  B of  Fig.  11a).  The 
maximum  impact  factor  remains  ousontially  the  same  in  both  cases. 
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Brldp;e3  of  More  Than  One  Floating  Support 

Linear  Behavior.-  Consider  the  tvro  ponton  bridge  sho\^(n  schema- 
tically in  Fig.  12. 


At  any  time  t the  unbalanced  force  at  any  one  of  the  concentrated 
masses  is  the  difference  between  the  applied  force  and  the  ponton  reaction 


^ 1-t  " ^1-t  " ^-t 


(28) 


IF  “ P ~ R 

^ 2-t  2-t  2-t 


For  a continuous  structure  of  this  type  the  reaction  of  either  ponton 
depends  on  the  displacement  of  both  pontons.  Assuming  linear  elastic 
behavioi-,  the  reaction  at  each  ponton  expressed  in  matrix  notation  is 


^12 

i^l-t 

^i-t 

^21 

(a22‘*^^) 

y2-t 

R2-t 

(29) 


Substituting  these  expi'cssions  for  the  reactive  forces  into 
equation  (28) 

^ 1-t  = ^1-t  -(^ii‘^*^)yi-t  - ^12y2-t 

^2-t  ^ ^2-t  ■ ^2i^l-t  "^^22^^^y2-t 


As  before,  I'rom  impulse  - momentum  considerations,  the 
integral  of  the  unbalanced  force  at  ponton  1 over  the  short 


time  interval  ut  is 
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^ Pj-l-  [(^ll^K)yi_i^a3^g72.i]  ■*-Pi_2-[a<  iyi-2~^  2~^^^y2-2] 


The  corresponding  change  in  momentum  of  ponton  1 is 
2ni2_ 


AV  = 


At  <^1-2  -^1-1  -''l-l 


Equating  impulse  to  change  in  momentum  the  expression  for  the 

displacement  of  ponton  1 in  terms  of  the  displacements  of  pontons 

1 and  2 at  time  t„  is 
4mi 


+^11  +K)yi_2  + ^12  72-2  " ^1 


where  the  load  coefficient  for  ponton  1 is 

Amp 

*Zl^ 


ifjui  ^ 

Ki  = P2^_;j_  +Pi_2  — TV — ^11  "'^)yi-l  "^12^2-1  .2  ^1-1^^^ 


Following  a similar  procedure  for  ponton  2 the  simultaneous 
linear  equations,  written  in  matrix  form,  are 


1 . Unu 

(-^+aii+K) 

^1-2 

^1-1 

&2i  ^“^■'^22^’^) 

^2-2 

“ 

^2-1 

or 


(30) 


|a'||  y]  ■ 1K| 

The  effect  of  virtual  iiiass  is  neglected  in  this  expression. 

To  express  the  displacements  y in  explicit  foria  it  is 
necessary  to  invert  the  matrix  la'|  . Letting  Itl  be  the 
inverse  matrix  equation  of  |a'j  equation  (30)  becomes 


T 11 

Ti2 

Kl-1 

yi-2 

y 21 

i^22 

^2-2 

or 


1^1 


•Ul 


(31) 
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A convenientpabular  scheme  can  then  be  used  to  succes- 
sively compute  the  displacements  of  the  pontons  at  any  time 
tjji  in  terns  of  the  displacements  and  velocities  at  time 
Fig.  13  shows  solutions  of  equations  (31)  different  times 

of  crossing.  The  physical  constants  for  the  two  ponton  bridge 
are  the  same  as  for  the  single  ponton  bridge  treated  earlier 
with  the  exception  that  the  sp^an  length  is  increased  to  1+5  feet. 

This  method  of  analysis  is  readily  extended  to  ponton 
bridges  with  more  than  two  pontons.  Fig.  14  shows  results  for 
a three  ponton  bridge  when  the  time  of  crossing  of  the  moving 
load  is  equal  to  Tq.  Virtual  mass  is  neglected. 

Non-Linear  Behavior.-  The  linear  problems  discussed  above  con- 
sist, essentially,  of  solving  the  linear  matrix  equation 
la'l  \y\  = U1 

which  would  allow  a direct  calculation  of  the  displacement  of 
each  ponton  frotn  the  matrix  equation 

\y\  = |K1 

If  virtual  nass  is  considered  the  solution  is  complicated 
by  the  fact  that  the  natrix  |a'(  is  displacement  decendeint. 

After  each  step  in  the  numerical  computation  the  diagonal  ele- 
ments of  the  matrix  |a'|  change  by  the  ai;iount 


4 


where  ^ is 


the  slope  of  tlie  virtual  iiiass-dispLacement 


-23- 


relationahip  at  time  and  is  the  difference  in  displace- 
ment over  the  time  interval  Jt,  = In  other  words,  |a'| 

is  changed  to 

la'I^  = laM  + Ibl  (33) 

where  \ bj  is  a diagonal  matrix. 

For  engineering  purposes  it  is  sufficient  to  consider  the 
change  in  |a'|  only  during  the  tirae  interval  when  cal- 
culations are  being  made  for  the  time  interval  This 

procedure  maintains  linearity  of  the  simultaneous  linear  equa- 
tions and  avoids  the  occurence  of  terms  of  an  order  higher  than 
one. 

The  problem,  therefore,  is  to  find  a new  inverse  matrix 

(34) 

It  would  be  laborious  to  invert  la'I'*''^  to  get  I 4^  1 ••  at 
each  step  in  the  computation.  Instead  it  is  exjjedient  to  express 
in  terii’.s  of  a ser'es  and  consider  only  the  first  few 
terms^”^^.  This  is  procedure  can  bo  used  only  if  the  corrections  Ibl 
are  small  relative  to  \a'|so  that  adequate  accuracy  is  obtained  from 
the  first  few  terms  in  the  series. 

From  equations  (33)  and  (34) 

- (la'l  + lbl)-l  (35) 

and  since 

= la'l  (l+la'l -^Ibl  ) 


(la'l  + Ibl) 


(36) 
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the:  ri^lit  hand  side  of  equation  (35)  becomes 
(i»'i  = {u'l  (1+ Ja'l -•‘•Ibl  )j 

- f(l!  la'I -Inn  (37) 

The  first  term  in  (37)  can  be  exi)ressed  as  the  continuous 
series 

^ [l+(-l)‘"  (»a'r^\bl)" 

n=l 

Therefore,  ^ 

HV/1^:-  = (la'I  + lbl)"^  = ^ [l+(-l)"(la'l"^lb»)'^]  la'l”^ 

n=l 

= Itl  - li^llblltH  +l4'nbll4'UbllYtv--- 

or 

h/|’^  =iirim^i  (38) 

where 

U 1 = I -Itllbl  (I'VlIbl  )2-r...  (39) 


Consider  an  exnjr.nie  of  this  procedure. 


4 

2 1 

0.3810 

-0.2860 

0.0476 

ja'l  = 

2 

11 

-0.2860 

0.7150 

-0.2860 

1 

2 4 

0.0476 

-0.2860 

0.3810 

4.1  2 1 

lO.l  0 

0 

2 3.1  2 

II 

0 

0.1 

0 

1 2 4.2 

0 

0 

r\  n 

9 

0.0381 

-0.0286 

0.0095 

-0.0286 

0.0715 

-0.0572 

0.0048 

-0.0286, 

0.0762 

"If  This  follows  fr*>iii  the  fundai 'cntal  rule  in  matrix  oi>eration 
that  (^^\  = \B\-^  U|-l 

L"  
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( Wllbl)^  = 


C.0023 

-0.C034 

0.0027 

C.0034 

C.0076 

-0.0087 

0.0014 

-0.0044 

0.0075 

0.964 

0.0252 

-0.0068 

0.0252 

0.9361 

0.0485 

-0.0034 

0.0242 

0.9313 

14^1*  = U'lKl  •= 


0.3598  -0.2558 
-0.2558  0.6482 
0.0361  -0.2480 


0.0361 

-0.2480 

0.3477 


The  correct  value  of  \4'(*  determined  by  inverting  |a'  1^ 


directly  is 


\r  I* 


0.3596 

-0.2552 

0.0359 


-0.2552  0.0359 
0.6467  -0.2472 
-0.2472  0.3473 


The  average  error  of  the  elements  of  the  iiatrixl^'/*  as 
computed  by  the  series  compared  to  their  exact  ■values  is  1.72^ 
vrith  a naximum  (at  t|/  22*)  of  2.20^. 

Considering  only  the  first  two  tenns  in  the  series 


0.9619 

0.0286 

-0.0095 

\yi  = i-ifiibi  = 

0.0286 

0.9285 

0.0572 

-0.0048 

0.0286 

0.9238 

0.3579 

-0.2519 

0.0340 

14'!’^  = WW^\  = 

-0.2519 

0.6393 

-0.2424 

0.0340 

-0.2424 

0.3436 
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The  average  error  considering  only  the  first  two  terms  of 
the  series  is  12.055^  with  a maximum  of  16.2^. 

Fig.  15  shows  the  results  of  applying  this  method  to  a 
three  ponton  bridge  for  different  times  of  crossing.  The  vir- 
tual mass,  for  each  change  in  the  )a'|  matrix,  was  determined 
for  each  ponton  from  Fig.  8.  In  general  it  is  not  always  necessary 
to  compute  a new  inverse  matrix  after  each  step  in  the  computa- 

tions. If  the  slope  of  the  virtual  mass-displacement  diagram  is 
small  it  may  only  be  necessary  to  correct  for  the  change  in  virtual 
mass  after  fcvery  second  or  third  step.  This  is  mostly  a matter  of 
Judgment  since  the  change  in  virtual  mass  (vMch  is  a function  of 
displacement)  during  one  step  in  the  computations  also  depends  on  the 
magnitude  of  the  moving  load.  The  corrected  inverse  matrix  for 

the  solutions  shown  in  Fig.  I5,  was  computed  after  every  second  step 
in  the  computations. 

Damping 

The  numerical  methods  described  above,  for  either  the  linear 
or  non-linear  problems,  are  easily  extended  to  include  the  effect 
of  viscous  damping.  The  damping  force  is  added  to  the  buoyant  resis- 
tance of  the  ponton  and;  from  equation  (12), 

IF  = l>(t)  - Ky  - CV  ikO) 

where  C is  the  coefficient  of  viscous  damping.  The  unbalanced  force 
over  a short  time  interval  becomes 

1 (t)  - Kyj^  - CV^  + l2(t)  - Ky^  - 


I 


2 
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Equating  this  impulse  to  the  change  in  momentum  and  rearranging 


terms 


(— ^+K)72  = ,\t2 


and  from  equation  (16) 

V.'ith  this  value  for  ^Q^^tion  (41)  becomes 


( 


_yi^+Kt-^)y  = (Jt«--K^-2^)y  + -^V  ^t  + P +P  (42) 

lilt‘d  at  ^ at2  at  at2  ■‘■2 


The  analysis  is  easily  extended  to  bridges  of  more  than  one 
ponton.  Fig.  16,  for  example,  showo  results  for  a three  ponton  bridge 

for  C°100  and  1000  Ib-sec/in.  The  logarithmic  decrement  corresponding 
to  these  values  of  C can  be  computed  from  the  eq\iation 

X = ili 

^cr 

where  ^ is  the  logarithmic  decrement 


'cr 


is  the  critical  damping 


The  critical  damping  can  be  computed  from 


„ ^ 4m  TT 

^cr  

To 

For  a throe  ponton  bridge neglecting  virtual  mass,  Tq  is  0.370  seconds 
and  Cj,j,  = 2060  Ib-sec/in.  The  logarithmic  decrements  corresponding  to 
C- " 100  and  1000  Ib-soc/in.  are  0,304  and  3.04  respect! velv. 


Accuracy  of  Numerical  1 rocedure 

The  stepwise  solution  of  the  initial  value  problems  discussed 
above  consists  of  the  evaluation  of  an  integral  at  pivotal  points 
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along  its  inter'/al  of  definition.  Consider  the  function  f(t)  shown 
in  Cig.'  17a.  The  integral  of  f(t)  from  t^^  to  tjjj+j  given  by 


’’f'd'l  . . fjfff^Ct)  + fjjj(t) 


f(t)dt  = 


jt  + e 


(43) 


m 


where  e is  the  small  shaded  area  in  Fig.  l?a.  For  a second  degree 
curve  this  error  is 


e = 2/^  lilt  (44) 

v/here  d is  the  maximum  vertical  ordinate  between  the  actual  curve  of 
f(t)  and  a straight  line  relationship  from  ^[^(t)  to 
derivation  of  the  numerical  procedure  used  above  it  was  assamed  that 
f(t)  varies  linearly  over  the  small  time  interval  at.  The  correc- 
tion e,  therefore,  indicates  the  discrepancy -between  the  correct  and 
the  approximate  values  of  the  solution  for  one  step  of  the  stepwise 
integration  process.  These  "truncation  errors"  accumulate  and  the 
results  after  a nximber  of  . teps  become  more  inaccurate  than  at  the 
origin. 


An  estimate  for  e,  however,  can  be  included  in  each  step  of 
the  numerical  integration.  If  this  estimate  is  conservative  enough 
to  include  the  largest  expected  value  of  e the  results  thus  obtained 
would  constitute  an  "upper  bound"  to  the  solution.  The  correct  solu- 
tion then  lies  between  the  values  obtained  by  including  e and  those 
obtained  by  neglecting  e. 

A safe  estijiate  for  the  ordinate  d between  the  actual  and 
assumed  curve  of  f(t)  over  the  time  interval  at  is  obtained  from 
Fig.  l?b.  Let  3 be  the  slope  of  f(t)  at  time  t . This  slope  can 

“I  lu 
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be  approximated  by 


Sm  = 


2^t 


The  di-tance  b of  Fig.  l?b  at  time  is  then 

_ o .eit  / _ ^ra+l(i)  “ 

b - Sm  /2 

4 


(45) 


(46) 


The  distance  c is 


c = 


^m+l(i)  “ ^m(i) 


(47) 


The  difference  between  b and  c is  a value,  d',  that  is  always 
larger  than  d for  a function  f(t)  of  any  degree.  The  error  e can 
then  be  approximated  by 

e'  =2/^  d*ot  (43) 

where  d'  = b - c 
From  equations  (46)  and  (47) 

<*'  * S - f.r.d*-)  - 


and 


= [2fm(t)  - W(t)  - Vl^^)j  ^^9) 


As  an  example  consider  the  single  i^ass  oscillator  shown  on 
page  9.  From  equation.;  (i?)  and  (14) 


f 


■'2 


(P(t)  - Ky)  dt  = muiV 


il 


or 


‘'2  (■‘^2 

r(t)dt  - K \ ydt  = m(V2-Vi) 


(50) 
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From  equation  (43) 


Vdt  = y^,  - y^ 


V2+V1 

= =-<at  + e. 


(51) 


and  from  equation  (49) 


= 


^ - V ''0]  (52) 

jubstituting  equation  (52)  into  equation  (51)  the  expression 


for  the  change  in  velocity  becomes 


h - h 


■(y^  - yp  - V2V1  + Vjv, 


^t  '•'2  •'1'  0 

The  integral  of  the  external  force  over  the  time  interval 
t2~tj^  becoiiies,  from  equation  (43  )> 


'2  Pl(t)  + P2(t) 

P(t)dt  >=  


(53) 


At  + ^ j^2P3^(t)  -P2(t)  -Po(t^ 

^ [2P2(t)  + 5Pi(t)  - Po(t)J  (54) 


and  in  a similar  manner 
't 


k(  ydt  = - + ^^(2y^  -y^  -y^) 


K£it 


r~  ^^2  ^5yi  -7q) 


(55) 


Substituting  equations  (53),  (54)  and  (55)  into  equation  (50) 
and  rearranging  tenr.s,  the  equation  for  the  displacement  of  a single 
mass  oscillator  at  time  t2>  in  t^rms  of  the  displacements  and  velo- 
cities at  tiii.es  tj^  and  tQ  and  the  applied  load^  becomes 
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(_l^+2K)y^  = (.i«m 
lX‘-  *- 


18m 


21m 


+ Kyo  + 2P2  + 5Pi  - Po 


(56) 


Numerical  Example.-  Again  considering  the  single  irass  oscillator 
shov/n  on  page  ||  and  assuming  that  m = 1 Ibs-sec  yin.  and  K = 

6/v  Ibs/in.,  the  fundamental  period  is 


0 


n 

T 


The  external  force  F is  chosen  to  vary  parabolical lly  with  time 
from  t=0  to  t=?Q,  or 

rU)  - 


4P'  .2 

~^-T  ^ 


‘0 


'0 


where  I ' = 1000  lbs. 

■'.'ith  these  numerical  values,  chosing  cit  = "^0/20>  equation 
(56)  becomes 


11823y2  = 11375/1  136A5V^E.t  - 19A9VQ<at 

+ + 2P2  " 5?!  - Fq  (57) 


The  displacement-time  relationship  for  the  single  mass  oscil- 
lator calculated  from  equation  (57)  is  shown  in  Fig.  18,  This 
solution  includes  the  corrections  e'.  When  the  corrections  a' 
are  neglected  equation  (19)  becomes 

2658/2  - 2530/1  ♦ 259/fViAt  + ?2  (58) 

Solutions  from  equation  (58)  are  shown  in  Fig.  18.  The  differ- 
ence between  the  maximum  ordinates  of  the  two  solutions  for  a 
time  interval  At  ■ ^0/20  is  1.13/^.  The  exact  solution,  shown 
in  Fig.  18,  lies  halfway  between  the  two  curves. 
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Fig.  18  also  shows  results  of  the  numerical  procedure  for 
^t  = "^O/lO*  The  difference  between  the  maximum  ordinates  for 
solutions  including  the  corrections  e^ and  for  solutions  neglecting 
this  correction  is 


ELECTRONIC  ATJALOG  SOLUTIONS 

The  solution  of  t*  i types  of  problems  discussed  above  by  nume- 
rical or  exact  means  is  laborious  and  time-consuming  when  a large 
number  of  problems  is  to  be  considered,  l^theiiatical  equations 
relating  measures  of  change  (derivatives)  of  this  ty]je,  however, 
are  particularly  amenable  to  solution  by  means  of  electronic  analog 
computors^^’^^  In  the:.e  computors,  the  variables  in  question  are 
represented  by  d-c  voltages  which  may  vary  with  time.  These  vol- 
tages are  made  to  obey  mathematical  relations  in  a manner  analogous 
to  those  of  the  original  problem.  Time  is  used  as  the  independent 
variable.  The  voltages,  or  machine  Vc  riables,  may  be  conveniently 
varied  and  measured  in  the  laboratory.  The  magnitudes  and  time 
scale  of  these  analogous  variables  can  be  adjusted  to  facilitate 
their  measurement.  The  solution  of  the  differential  equations  of 
motion  of  the  pontons  proceeds  as  follows: 

a)  VJhen  the  machine  is  set  up  to  solve  the  problem  the 
voltages  (machine  variables)  are  set  to  the  initial 
conditions  of  the  given  problem. 

b)  The  coraputor  is  then  started  and  forces  the  voltages 
to  vary  in  the  manner  prescribed  by  the  differential 
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II 


equations.  The  output  voltaj^es  are  ^.lotted  as  a func- 
tion of  tijne  and  constltutt  a solution  of  the  {problem. 

A more  detailed  description  of  the  basic  con;ponents  and  pro- 
gramming is  given  in  Appendix  C, 

Solutions  were  obtained  by  the  analog  comjjutor  for  the  following 
cases: 

A.  Single  ponton  bridge  con;. i deling  the  i.:ass  of  the  moving 
load  P.  The  deflected  load  line  is  approxint^ted  by  a second 
degree  curve  through  the  sup>)orts  and  the  single  ponton. 

(1)  Neglcct.ing  virtual  nass  an<i  damping,  = 3m  and  qm 

(2)  Considering  virtual  mass,  n.j^  = 3m 

B.  Two  ponton  bridge  neglecting  the  mass  of  the  moving  load  P, 

(1)  Neglecting  virtual  nass  and  dam:ping. 

(2)  Considering  virtual  mass. 

(3)  Considering  virtual  nass  and  damping,  .solutions  were 
obtaimd  f f r d:>.i;f>ing  coefficients  of  100  and  lOOC  .Ib-sec/in. 

C.  Three  j onton  bridge  neglecting  the  i £iss  of  the  moving  load  P. 

(1)  I'egl  ccti'ig  virtual  mass  and  dar.iplng  .and  choosing  the 
value  of  EI/l3  - R5.8  Ibs/in. 

(2)  The  san-.e  as  (l)  excep't  EI/L^  = A2.9  Ibs/in, 

(3)  Considering  virtual  mass  and  damping.  The  d.amping 
coefficients  were  assumed  to  be  100  .and  1000  Ib-sec/in. 

D.  Five  i onton  bridge  neglecting  the  nass  of  the  n.ov  Lng  load  P. 

(1)  iJeglecting  virtual  i.aos  and  d-aiiiping. 

In  all  the  case:  described  above  soluti  ns  were  obtained  for 
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different  tines  of  crossing. 

Results  of  Analog  Solutions 

Case  A-1.-  Since  the  shape  of  the  deflected  load  line  is  assumed 
the  displacement  of  the  bridge  at  any  jioint  x is  known  in  terms 
of  the  displacement  of  the  p-onton 

and  the  vertical  acceleration  of  the  bridge  at  point  x,  neglecting 
the  Coriolis  acceleration,  becomes 

Vx  = e(x)yp  (59) 

Since  the  speed  of  the  ireving  load  is  constant  the  position 
of  the  load  can  be  exiressed  in  terms  of  time  t and  equation  (59) 
becomes 

= 9(t)yj,  (60) 

For  a second  degree  deflection  curve 


e(t) 


(61) 


The  force  of  the  moving  load  i ap>plied  to  the  bridge  at  the 
point  X is 


F 


X 


(62) 


substituting  equation  (60)  into  equation  (62) 

^x  = ^ ® ^p 

The  reaction  -'^(t,yp)  on  the  ponton  due  tc  the  force  is 
assujiied  to  vary  linearly  whence, 

= F^0(t)  = [p-m^  e(t)Vp]  0(t)  (6A) 


F 
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Th&  function  0(t)  represents  the  influence  of  the  force 
Fx  at,  any  time  t on  the  ponton.  If  H is  the  time  for  the  load 
to  cross  the  bridge  this  function  can  be  represented  as  shovm 
in  Fig.  19. 


Using  equations  (6l)  and  (64)  the  equation  of  motion  for 
the  ponton  of  a single  jonton  bridge  (n=l),  conrddering  the  nass 
of  the  load  and  neglecting  virtual  mass  and  dai.ijdng,  becomes 
from  equation  (4) 

Introducing  the  follov/ing  numerical  vulues 
= 3m  = 60.6  Ibs-sec^/in. 
aii+K  = 34422  Ibs/in. 

P = 23,400  lbs 

El  = 320  X 10®  lbs-in.2 

L = 30  ft 
Tq  = 0.152  sec 

equation  (65)  becomes 

y + 17047  = [1153.4  + (^1^ yj  0(t)  (66) 

Solutions  for  ciouation  (66),  taken  fro'r  the  computor  record, 
••'.re  -ihi  .-/n  in  Fig.  20  for  H=0.5,  1.0,  1,5  2.0  Tq, 
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Aiisuirdr.^-;  nij^  = 5n-  tcucition  (65)  becomes 

j ♦ 1704y  4i931  + ^)y^0(t)  (67) 

L H J 

Computer  solutions  for  equation  (67)  are  shovm  in  Fig.  21 
for  H = 0.5,  1.0,  1..5,  2.0  and  2.5  Tq. 

Case  A-2.-  From  Fig.  8 the  virtual  mass,  at  the  dead  weight  dis- 
placement of  the  ponton  (5.2)inches),  is 

m'  = M + J (y)  = 74  + 2.23y  (68) 

With  this  value  and  using  eiipression  (64)  equation  (6), 
neglecting  damping,  becomes 

(74^2.23y)y  + 34422y  = [_23400-60.6(-^i (^9) 

Solutions  for  equation  (69)  are  given  in  Fig.  22  for 
H =0.5,  1.0,  1.5,  2.0  and  2.5  Tq. 

Case  B-1.-  From  equation  (2)  the  equations  of  motion  of  a two 
ponton  bridge  (n=2),  with  the  following  numerical  constants: 

P = 23,400  lbs  L = 45  ft 

El  = 320  X 108  lb-in. 2 
K = 1500  Ib/in. 
m = 20.2  lbs-sec2/in. 
becomes 

20. + 54l75yi  - 46O90yp  = F.(t)  = 23400  0,(t) 

■L  J-  ^ X 1 

20.272  - 460907^^  + 5417572  = ^^(t)  - 23400  0^(t) 

The  functions  0j^(t)  and  02(t)  represent  the  influence  of 
the  moving  load  1',  at  any  time  t,  on  pontons  1 and  2,  respectively 
These  functions  can  be-  represented  as  shown  below  in  Fig.  23. 
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FIG.  23  INFLUENCE  OF  MOVING  LOAD  ON  FONTONS  1 AND  2 
Solutions  for  equations  (70)  are  given  in  E’ig.  24  for 

different  times  of  crossing. 

Case  B-2.-  V/ith  the  value  for  the  virtual  mass  given  in  expres- 
sion (68)  numerical  values  of  Case  B-1  equations  (3)  become 

(74  + 2.23yi)y,  + 54l75y.  - 46C90y„  = F.  (t)  = 23400  0, (t) 
1X21  1 

(74  + 2.23y2)y2  ~ 54l?5y2  = 

Fig.  25  shows  the  computor  solutions  for  equations  (71) 

for  different  tir  es  of  cros  ing. 

Case  B-3.-  Considering  virtual  vass  and  damping  and  neglecting 
the  mass  of  the  moving  load  equations  (6),  with  the  numerical 
values  of  Case  B-1,  become 

(74  + 2.23y^)>\  + Cy^  ->  54175y3_  - = F^(t)  = 23400  0j_(t) 

( 

(74  + 2.23y2)>2  ^ + 54175y2  = F^(t)  = 23400  0^(t) 


Figs.  26  and  27  show  solutions  for  equations  (?2)  for  C =100 
and  1000  Ib-sec/in.,  respectively. 

Case  C-1.-  For  a three  ponton  bridge  (n=3),  neglecting  virtual  mass 
and  damping  and  using  the  same  numerical  values  of  Case  B-1  with  the 
exception  that  L=60  ft,  equations  (2)  are 
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20.2^^  + 55000yj^  - 493007^  + 2O6OCV3  = 23400  0^(t) 

20.272  - 4930073^  + 7550072  “ ^93007^  = 23400  0^(t)  (73) 

20.273  + 206007^  - 493OO72  ■+  550007^  = 23400  0^(t) 

The  functions  0 (t)j  0p(t)  and  0, (t),  \rtiich  reprocent  the 

effect  of  a unit  load  at  anv  time  t on  pontons  1,  2 and  3> 
respectivcl7,  can  be  represented  as  shovm  in  Fig.  28. 


Fig.  29  gives  solutions  for  equations  (?3)  for  different 
times  of  crosring. 

Case  C-2.-  This  case  is  the  same  as  Case  C-1  with  the  exception 

O p 

that  the  vjilue  of  El  is  reduced  from  320  x 10  Ibs-in.'^  to 
160  X 10®  Ibs-in.^.  ’With  this  change  equations  (2)  become 
20.^1  + 282507^  - 2465072  + 103007^  = 23400  0^(t) 

20.272  - 2465073_  + 385OO72  - 2465O73  = 23400  02(t)  (74) 

20.^3  + 103007^  - 246507^  + 282507^  = 23400  03(b) 

Fig.  30  gives  the  computer  solutions  for  equations  (74) 
for  a range  of  crossing  timeu. 

Case  C-3.-  Considering  virtual  trass  and  damping  and  the  nume- 


rical values  from  Case  B-1  equations  (6),  neglecting  the  mass 
of  the  moving  load,  become 
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{7h-^2.23y^)\  + + 5500Qyj^  - 493007^  + 20600y^  = 23400  0^{t) 

(74+2.2372)72  + Cjr^  - 493007^  + 7^^00y^  - 493007^  = 23400  02(t)  (75) 

(74+2.2373)2^3  + 07^  + 206007^  - 493007^  + 55OOO73  = 23400  0^(t) 

Figs.  31  aJid.  32  present  solutions  for  equations  (75)  for 
C=100  and  1000  Ib-soc/in.,  respectively,  and  for  times  of  cross- 
ing of  3.0,  4.0  and  6.0  T^. 

Case  r— 1.-  Using  the  numerical  values  from  Case  B-1  equations 
(2)  for  a five  ponton  bridge  (n=5),  become 
20. 2^^  + 55750y^  - 52240y^  + 223h5y^  - 6114  7^  +1540y^  = 23400  0^^(t) 

20.^2  - 522407^  + 7860072  - 583407^  + 24370y^  -611475  = 23400  02(t) 

20.^3  + 2284571  - 583hOy^  + 8013073  - 583407^+2284575  = 23400  0^(t)(76) 
20. 2y^  - 6114  7i  + 243707^  - 6834073  + 78600y^-52240y^  = 23400  0^(t) 

20.^*^  + 1540  7i  - 6114  72  + 2284573  - 52240y^+ 5575073  = 23400  05(b) 

Fig.  33  shows  computer  solutions  of  equations  (76)  for 
different  times  of  crossing. 

Accuracy  of  Analog,  Solutions 

The  electronic  analog  coinputor?;  described  above  are  essentially 
electrical  devices.,  '.hen  these  dev.ic  cs  are  "set  uj' " to  solve  a set 
of  di  ffercntial  equations  they  :irt-  wired  in  such  a way  that  the 
e(ii.itionr.  describ.'tng  tiic  variati.  n of  voltage  with  time  are  analo- 
gous to  the  set  of  Ji  f fc^renti al  equations  \/hose  solution  is  d'Vvired. 

The  con.:tantr  of  these  equations  .are  represented  in  the  computor  by 
resistances  wiri  ch  are  .odjiisted  to  tht:i  r correct  value  by  rheostats. 

The  .'accuracy  of  the  cf.m;utcr  (t  the  degree  to  whi  on  n record  of  vol- 
tage van.',  ation  in  the  coii:putor  can  be  said  to  r'ciu'esent  a solution 
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of  the  differential  equations)  is  determined  by  the  accuracy  of 
these  rheostats  and  the  reliability  of  its  components.  Tl.ese  com- 
ponents perform  the  necessary  matheiratical  operations  of  addition, 
multiplication,  division,  and  integration. 

As  a check  on  this  accuracy  the  exact  solution  of  equation  (70) 
describing  the  behavior  of  a two  ponton  bridge  was  comi>ared  to  the 
analog  solution.  These  results,  for  ponton  1,  are  shewn  in  Fig.  34- 
There  is  no  difference  in  the  maximum  value  of  displacement  in  each 
case.  After  the  peak,  however,  the  analog  solution  lags  the  exact 
solution  by  a small  amount. 

The  numerical  solution,  from  equation  (31)  is  also  shown. 

1I'!FULSB-DIS:  LACIMENT  TESTS 

Experimental  imj'ulse-displacenient  tests  were  conducted  on  pontons 
of  different  shapes  in  an  open  58-in.  flujne.  The  pontons  were  sub- 
jected to  an  arbitrary  impulse  which  was  measured  and  recorded.  The 
corresponding  displacements  were  also  measured.  These  displacements 
v/ere  then  compared  to  calculated  displacements  based  on  the  known 
impulse  to  afford  a measure  cf  the  reliability  of  the  virtual  mass 
and  damping  factors.  Several  assumptions  are  made  in  computing  the 
displacements  from  the  given  impulse  by  the  numerical  procedure  pre- 
viously discussed.  It  is  assumed  that  the  water  level  does  not  change 
during  the  displacement  - that  is  to  say,  that  the  transient  waves 
formed  at  the  sides  of  the  pontons  due  to  the  rapid  displacement  of 
water  are  neglected.  It  is  also  assumed  that  the  buoyant  force- 
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disi'lace.’iient  relationship  is  linear  and  equal  to  the  slope  of  the 
buoyant  force-displacement  curve  at  the  jioint  of  static  dead  weight 
displacement.  The  virtual  mass  at  each  step  in  the  computations  is 
determined  from  the  experimental  virtual  ma-s- displacement  curve  for 
each  ponton. 

The  experiiuental  technique  for  determining  the  added  masses  for 
the  models  is  based  on  a method  developed  by  T.  E.  5telson^^*^^ . This 
method  consists  of  connecting  the  object  to  be  tested  to  a metal  rod 
which  is  rigidly  co'.inected  to  the  center  of  a simple  beaun.  The  beam 
is  then  vibrated  with  the  object  in  the  air  and  again  when  it  is 
submerged  to  the  desired  depth  of  water. 

The  frequency  of  a simple  beam  with  a weight  concentrated  at  the 
center  is  given  by 


v/here  f = frequency,  (cycles/sec) 

K = spring  constant,  (ibs/in.) 

m = mass  of  suspended  object  and  equivalent  concentrated 

weight  of  beam  (one-half  total  nass  of  beam),(lbs-sec2/in. ) 

The  ratio  of  the  trass  of  a ponton  in  water  to  its  mass  in  air 

can  then  be  calculated  from  the  relationship 

_ fa^ 
fw^ 

where  m^  = apparent  mass  of  ponton  submerged  in  water 
= mass  of  ponton  in  air 

fa  = frequency  of  bean,  and  pcnton  when  vibrated  in  air 
fw  = frequency  of  beam  and  ponton  when  vibrated  with  ponton 
at  desired  subniergence 


The  test  beam  is  desigaed  to  have  a relatively  hi^  frequency 
(20-30  cycles/sec)  within  the  expected  range  of  virtual  masses.  By 
keeping  the  amplitudes  small  the  ratio  of  inertia  forces  to  viscous 
forces  is  high^^^\ 

The  apparatus  used  to  determine  the  added  masses  is  shown,  in 
Fig.  35  and  schematically  in  Fig.  36.  The  test  beam  is  suspended 
longitudinally  between  two  I-beams  which,  in  turn,  are  cast  in  heavy 
concrete  blocks.  This  supporting  framework  has  a high  frequency 
with  respect  to  the  natural  frequency  of  the  best  beam  to  prevent 
the  occurrence  of  harmonics  or  resonance  which  would  distort  the  test 
results.  A magnetized  rod  is  fixed  vertically  near  the  center  of 
the  test  beam.  A detector  coil  is  suspended  around  this  rod  and 
connected  to  a Brush  Oscillograph.  As  the  test  beam  vibrates  an 
alternating  emf  is  set  up  in  the  coil  and  is  recorded  on  the  oscillo- 
graph. The  frequency  of  vibration  is  determined  from  the  oscillo- 
graph chart. 

The  three  types  of  model  pontons  tested  are  shown  in  Fig.  37 
and  schematically  in  Fig.  38.  The  corresponding  bouyant  force-displace- 
ment relationships  are  shown  in  Fig.  39.  In  order  to  facilitate  the 
testing  procedure  the  test  beam  is  first  calibrated  - that  is,  known 
weights  are  added  to  the  beam  and  the  natural  frequency,  corresponding 
to  each  weight,  is  noted.  The  curve  thus  obtained  is  shown  in  Fig.  10. 
The  ponton  to  be  tested  is  then  fixed  to  the  beam  and  the  corresponding 
natural  frequency  noted  on  the  curve,  btien  the  ponton  is  then  vibrated 
in  water  the  effective  added  weight,  due  to  the  acceleration  of  the 
ponton  through  the  water,  is  given  by  the  difference  between  the 
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the  points  on  the  abscissa  corresponding  to  the  two  different  fre- 
quencies. The  results  of  these  tests  are  shown  in  Fig.  41  which 
gives  the  virtual  uiass  for  the  model  pontons  as  a function  of  dis- 
placement. 

A similar  procedure  was  followed  for  determining  the  virtual 
mass-displacement  relationship  of  an  M-4  ponton.  A one-twelfth 
scale  model,  shown  in  Fig.  42,  was  used  for  this  purpose.  The  tests 
were  carried  out,  in  this  case,  in  an  18-in.  flume. The  added  mass, 
at  any  particular  displacement,  is  expressed  in  terms  of  the  iiass  of 
displaced  water.  The  added  mass  of  the  full  scale  ponton,  at  the 
corresponding  full  scale  displacement,  is  determined  by  multiplying 
the  added  mass  of  the  model  at  that  displacement  by  the  ratio  of  the 
weight  of  displaced  water  of  the  full  scale  ponton  to  the  weight  of 
displaced  water  of  the  model  ponton.  By  this  method  the  virtual 
mass-displacement  relationship  shown  in  Fig.  8 was  determined. 

The  aj)jjaratus  used  for  the  impulse-displacement  tests  is  shown 
in  Fig.  43  and  schematically  in  Fig.  44.  The  vertical  rod  attached  . 
to  the  center  of  the  ponton  is  restricted  to  vertical  motion  by  the 
slide  bearing  atta  ,hed  to  the  side  of  one  of  the  I-beams  used  in 
the  virtual  nass  apjaratus . The  force  on  the  ponton  is  determined 
from  SR-4  electric  strain  gages  attached  to  an  aluminum  "horseshoe" 
at  the  center  of  the  j’onton.  The  force  is  applied  to  the  vertidal 

^'-The  proximity  of  the  slvie  walls  has  an  effect  on  the  virtual 
mass  ccef fJ cient.  The  values  used  for  a study  of  a full  scale  j;onton, 
however,  are  t .nly  lucmt  Ij  be-  re>rerentatJ  ve  and  no  effoi  t was  ii.ade 
tc  cr.Tvct  fv-rthib. 
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rod  which  is  attached  to  the  top  o£  the  "horseshoe"  and  bending; 
stresses  are  induced.  The  "horseshoe"  is  calibrated  so  that  the 
vertical,  force  in  the  bar  can  be  determined  from  the  strain  in  the 
lR-4  c^ces  due  to  these  bending  stresses. 

Displacements  are  determined  in  much  the  same  manner  by  a 
cantilevered  strip  of  thin  spring  steel  to  which  SR-4  strain  gages 
are  fastened.  The  sp-ring  is  calibrated  so  that  the  vertical  dis- 
placement of  the  free  end  can  be  interpreted  in  ter.ms  of  the  strains 
in  the  SR-4  gages.  The  end  of  the  steel  strip  is  deflected  by  a 
metal  "finger"  attached  to  the  ponton.  The  spring  is  thin  enough 
so  that  its  resistance  to  displacement  of  the  ponton  is  negligible. 

Fig.  45  shows  results  of  an  iinpulse-disi)lacement  test  on  the 
rectangular  ponton.  The  dLS]lacciTent  curves  .calculated  from  the 
recorded  impulse  with  virtual  .rass  considered  and  with  virtual  mass 
and  damping  considered  agree  well  i«.th  the  experin.ental  displacement 
curve.’!-  The  calcul-  ted  di.s;  lacer.ent  curve  with  the  virtual  mass 
neglected,  however,  shows  practically  no  agreement  with  the  experi- 
mental results.  An  inspection  of  Fig.  45  indicates  that  in  this  case 
the  -virtual  mass  is  of  greater  impiortance  for  design  considerations 
than  damping.  Neglecting  damping  gives  results  on  the  safe  side 
wl-iereas  neglecting  -virtual-  mass  yields  results  considerably  on  the 
unsafe  side.  Siiiiilar  good  agreement  between  calculated  and  experi- 
mental displacemient  curves  for  the  rectangular  ponton  is  shown  in  Fig,  46. 

!!The  determination  of  the  damping  coefficient  is  shown  in 


Appendix  E. 
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Fig.  47,  48,  and  49  present  results  of  impulse-dis placement 
tests  on  the  ponton  with  triangular  cross-section.  Figs.  47  and  48 
indicate  the  importance  of  damping  on  tests  of  short  dxaration  with 
relatively  high  velocity.  The  relatively  large  error  (I'i^C)  between 
the  calculated  displacement  curve  with  damping  considered  and  the 
experimental  displacement  curve  i i Fig.  48  can  probably  be  attributed 
to  inaccuracies  in  .etermining  the  damping  coefficients. 

Results  of  tests  on  the  ponton  with  semi-circular  cross-section 
are  shown  in  Figs.  50,  51  and  52.  The  agreement  between  calculated 
and  experimental  displacement  curves,  for  the  test  of  short  duration 
shown  in  Fig.  52,  is  not  as  good  as  in  the  case  of  the  other  pontons 
tested. 

SUIOIARY 

The  relationship  between  maximum  impact  factor  and  number  of 
pontons  is  shown  in  Fig.  53  Tor  bridges  of  one,  two,  three  and  five 
pontons.  The  open  circles  shou  results  neglecting  virtual  mass, 
damping,  and  the  mass  of  the  load.  Tiio  scatter  in  values  of  impact 
factor  vvith  increasing  n\jmber  of  pontons  is  not  unreasonable.  Consider 
the  mici-span  cieflection  of  identical  simple  span  beams  of  uniform 
section  subjected  to  different  loads  which  produce  the  same  maximum 
fiber  stress  at  mid-span.  The  ratio  of  the  mid-span  deflection  of 
a beam  subjected  to  equal  end  moments  to  the  deflection  of  a beam 
under  uniform  load  is  1.20.  Similar  ratios  for  beams  with  2 and  4 
equal,  symmetrically  placed  concentrated  loads  is  1.02  and  l.Ol.  For 
beams  with  1,  3 and  5 equal  symmetrically  placed  concentrated  loads 


the  ratios  are  0,80,  0.95  and  0.98.  As  the  number  of  load  concentra- 
tions increase  the  ratio  approaches  the  limiting  value  of  1,0.  The 
ratio^  for  odd  number  load  concentrations,  approaches  the  limiting 
value  asymptotically  from  below;  while,  for  even  number  load  concen- 
trations, the  approach  is  from  above.  The  plotted  data  of  Fig. 53 
indicate  this  trend.  The  variations  may  be  attributed  to  the  contri- 
bution of  the  higher  modes  of  vibration  to  the  displacements  in  each 
case. 

Fig.  54  shows  the  first  and  second  modes  of  vibration  for  the  four 
bridges  summarized  in  Fig.  53.  For  a one  ponton  bridge  only  the  first 
mode  is  excited  by  the  action  of  the  moving  load  P.  Fib,  54b  shows 
hov/  the  first  two  modes  are  excited  in  a two  ponton  bridge  by  the 
moving  load  and,  moreover,  how  all  of  the  concentrated  mass  of  the 
bridge  oscillates  in  this  second  mode.  One  Would  then  expect  the  dy- 
namic displacements  to  be  greater  for  a two  ponton  bridge  than  a one 
ponton  bridge  due  to  the  contribution  of  the  second  mode.  This  con- 
tention in  supported  by  Fig,  53. 

In  the  case  of  a three  ponton  bridge,  however,  only  66,6'^  of  the 
concentrated  nvass  oscillates  in  the  second  mode.  The  center  mass  is 
at  a node.  One  v/ould  then  expect,  neglecting  the  contribution  of  the 
third  mode,  that  the  rmiximum  dynamic  displacements  of  a three  ponton 
bridge  would  be  less  than  a two  ponton  bridge  but  more  than  a one 
ponton  bridge.  This  is  supported  by  Fig.  53.  For  a five  ponton  bridge 
8C^  of  the  concentrated  mass  oscillates  in  the  second  mode  as  is  shown 
in  Fig.  54  d and  the  maximum  dynamic  displacements  are  greater  t’nan 
those  of  a tnree  ponton  bridge  (66.6?^  of  the  oscillating  mass  in  the 
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second  mode)  and  less  than  those  of  a two  ponton  brid/^e  (100  of 
the  mass  oscillating  in  the  second  mode). 

As  the  number  of  pontons  increases  the  difference  in  the  impact 
factor  between  any  two  consecutive  values  of  n will  decrease  and 
eventually  the  relationship  between  impact  factor  and  n will  be 
consistent. 

Fig.  10  shows  the  impact  factor  as  a function  of  time  of  crossing 
for  a single  span  bridge  considering  virtual  mass  and  neglecting 
virtual  mass.  The  maximum  ordinate  is  the  same  in  each  case.  The 
effect  of  the  virtual  mass  is  to  cause  the  maximum  impact  factor 
to  occur  at  a slower  speed  of  crossing.  Figs.  25a  and  29a  show 
similar  comparisons  for  two  and  three  ponton  bridges,  respectively. 

The  maximum  impact  factor  for  a two  ponton  bridge  is  increased  when 
virtual  mass  is  considered.  In  the  case  of  the  three  ponton  bridge 
a consideration  of  virtual  mass  leads  to  lower  impact  factors. 

These  results  are  also  summarized  in  Fig.  53.  The  open  squares 
are  results  for  bridges  of  one,  two  and  three  pontons  considering 
virtual  mass. 

The  effect  of  the  mass  of  the  moving  load  was  considered  for 
a single  ponton  bridge  only.  The  impact  factor-time  of  crossing 
relationship  shown  in  Fig,  11  for  different  values  of  the  ratio 
mj^m  indicates  that  the  maximum  impact  factor  is  the  same  for  all 
values  of  this  ratio.  The  deflected  load  line,  in  this  case,  was 
assumed  linear.  Fig.  55  shows  similar  results  for  mL  * 3m  and  5m 
which  are  summarized  in  rigs,  20  and  21.  In  this  case,  however,  the 
deflected  load  line  was  assumed  to  be  a second  degree  curve.  The 
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maxlmum  impact  factors,  however,  are  the  same  as  those  of  Fig.  11 
indicating  that  the  assumption  of  a linear  deflection  curve  v lelds 
sufficiently  accurate  results  when  the  mass  of  the  load  is  taken 
into  account.  This  is  emphasized  by  Fig.  55  which  shows  a computor 
solution  for  a single  ponton  bridge  consideidng  virtual  mass  and  the 
mass  of  the  moving  load  assuming  a second  degree  deflection  curve 
(Fig.  22)  and  a numerical  solution  of  the  same  problem  assuming  a 
linear  deflection  curve.  Fig.  55  also  indicates  that  the  considera- 
tion of  virtual  mass  leads  to  a lower  impact  factor  for  a single 
ponton  bridge,  when  considering  the  mass  of  the  lead,  than  is  obtain- 
ed when  virtual  mass  is  neglected. 

The  effect  of  varying  the  ratio  of  K to  EI/l3  for  a three  ponton 
bridge  is  shown  in  Fig.  57  which  is  a plot  of  impact  factor  as  a 
function  of  KL^/EI.  The  results  for  values  of  KL3/EI  of  17.5  and 
35  are  obtained  from  Figs.  29  and  30*  A value  for  KL3/ei  of  infinity 
means  that  there  is  no  elastic  interaction  between  the  pontons.  This 
implies  that  each  ponton  and  its  adjacent  spans  act  as  single  ponton 
bridges.  The  impact  factor,  therefore,  from  Fig.  53>  is  1.52. 

A value  of  KL3/EI  of  zero  implies  that  the  value  for  K is  zero. 
The  equations  for  a three  ponton  bridge,  for  K“0,  are  then  the  same 
as  equations  (lO)  with  n“3«  Solutions  of  these  equations  for  a 
moving  load  P and  for  different  times  cf  crossing  are  given  by  J.  M. 
Bonetti.^^^  The  maximum  impact  factor  obtained  is  1.67. 

The  effect  of  damping  on  impact  factor  ip  shown  in  Fig.  58  for 
a tvro  ponton  bridge  and  Fig.  59  for  a three  ponton  bridge.  These 
plots  show  the  relationship  between  impact  factor  and  time  of  crossing 
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for  damping  coefficients  C ■■  0,  100  and  1000  Ib-sec/in,  Th*a«  curves 
are  summarized  in  Fig.  60.  Curve  C of  Fig.  59  reveals  that  the  im- 
pact factor  approaches  1.0  at  large  times  of  crossing  for  C^IGOO 
Ib-sec/in.  Curve  B of  Fig.  60,  which  shows  the  relationship  between 
impact  factor  and  coefficient  of  viscous  damping,  for  the  three  ponton 
bridge,  then  has  the  value  of  1.0  at  C-  1000  Ib-sec/in.  and  also  has 
zero  slope  since  all  the  values  of  C above  1000  Ib-sec/in.  (and  some 
below)  cause  an  impact  factor  of  1.0. 

CONCLUSIONS 

1.  The  numerical  procedure  discussed  above  can  be  satisfactorily 
used  for  the  analysis  of  the  dynp.mic  response  of  floating 
bridges  to  transient  loads.  The  accuracy  obtained  compares 
well  with  that  obtained  by  an  electronic  analog  computor  and 
is  sufficient  for  engineering  purposes.  The  method  is  readily 
extended  to  include  damping,  the  mass  of  tho  moving  load  and 
non-linear  behavior. 

2.  Virtual  mass  has  a significant  effect  on  the  response  of  floating 
bridges  to  transient  loads.  The  effect  is  more  pronounced  than 
the  effect  of  viscous  damping.  Moreover,  the  omission  of 
virtual  mass  in  the  analysis  of  a floating  bridge  leads  to 
displacements  on  the  unsafe  side. 

3.  The  impact  factor  for  a single  ponton  bridge  does  not  change  when 
the  mass  of  the  load  is  taken  into  account.  Varying  the  ratio 

of  the  mass  of  the  moving  load  to  the  mass  of  the  bridge  merely 


causes  the  maximuin  displacement  to  occur  at  a different  speed 
of  crossing. 

The  assumption  of  a linear  deflection  curve,  v«fhen  considering 
the  mass  of  the  load,  is  a good  approximation  for  the  shape  of 
the  deflected  structure.  Comparison  of  solutions  for  a single 
ponton  bridge  considering  a linear  deflection  curve  and  consi- 
dering a second  degree  (parabolic)  deflection  cvu*ve  show  complete 
agreement . 

The  ratio  of  the  buoyant  resistance  K of  the  ponton  to  the  stiff- 
ness EI/l3  of  the  bridge  has  some  effect  on  the  impact  factor. 

For  higher  values  of  this  ratio  the  impact  factor  decreases 
and,  neglecting  virtual  mass  and  the  mass  of  the  moving  load, 
approaches  an  impact  factor  of  1.52.  The  upper  limit  of  the 
impact  factor  for  lower  values  of  this  ratio  is  1.6?. 

The  virtual  mass  and  damping  coefficients  obtained  by  experi- 
ments are  reliable.  Calculated  displacements,  based  on  these 
values,  compare  well  with  experimental  results. 
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SINGLE  rONTCN  BRIDGE 


Eauation  of  f-^otion: 


(^,K)y2  . (-^-K)yi  . 


H=Tq 

^t=To/20 

1=23,400 

K=34,420 

Time 

Pl^P2 

To 

y 

Vdt 

0 

0 

0 

0 

0.05 

2,340 

0,00164 

0.00328 

0.10 

7,020 

0.00967 

0.01278 

0.15 

11,700 

0.02986 

0.02760 

0.20 

16,330 

0.06682 

0.04632 

0.25 

21,060 

0.12354 

0.06712 

0.30 

25,740 

0.20110 

0.08800 

0.35 

30,420 

0.29858 

0.10696 

o.z*o 

35,100 

0.41314 

0.12216 

0.Z*5 

39,780 

0.54029 

0.13214 

0.50 

14,,U60 

0.67433 

0.13594 

0.55 

44,460 

0.80562 

0.12664 

0.60 

39,780 

0.91824 

0.09860 

0.65 

35,100 

0.99480 

0.05452 

0.70 

30,420 

1.02138 

-0.00036 

0.75 

25,740 

0.98985 

-0.06270 

0.80 

21,060 

0.89572 

-•  .12556 

0.85 

16,380 

0.74150 

-C. 18288 

0.90 

11,700 

0.53550 

-0.22912 

0.95 

V ,020 

0.29102 

-0.25984 

1.00 

2,340 

0.02506 

-0.27208 
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single  lONTON  BRIDGE  CX3NS1DERING  VIRTUAL  NASS 


AH  LIEU  LOAD-TIMl  CURVE  FOR  R-T^ 

Equation  of  Motion; 


m ^ 

(40. 568-up 

-i-T 'll.  = ( 

'40  568 

m , 

K' 

or 

Ay2  = + 

cvi^^t  t : 

D 

m'  = 74.-''2.23y 

c.t 

= ^®/20 

Time 

^0 

A 

0 

c 

D 

y 

V^t 

0 

149.613 

147.613 

148.613 

0 

0 

0 

0.05 

149.613 

147.613 

148.613 

0.67984 

0.00454 

O.OO9O8 

0.10 

149.633 

147.633 

148.633 

2.03953 

0.02713 

0.03610 

0.15 

149.734 

147.734 

148.734 

3.39920 

0.08533 

O.O8O3O 

0.20 

149.995 

147.995 

148.995 

4.75888 

0.19568 

0.14040 

0.25 

150.489 

148.489 

149.489 

6.11856 

0.37320 

0.21464 

0.30 

151.284 

149.284 

150.284 

7.47824 

0.63092 

0.30080 

0.35 

152.438 

150.438 

151.438 

8.83792 

0.97945 

0.39626 

O.AO 

153.999 

151.999 

152.999 

10.19760 

1.42664 

0.49812 

0.45 

156.002 

154.002 

155.002 

11.55728 

1.97736 

0.60332 

0.50 

158.468 

156.468 

157.468 

12.91696 

2.63343 

0.70882 

0.55 

161. 406 

159.406 

160.406 

12.91696 

3.38'^26 

0.79484 

0.60 

164.773 

162.773 

163.773 

11.55728 

4.20433 

0.84330 

0.65 

168.441 

166.441 

167.441 

10.19760 

5.05324 

0.85452 

0.70 

172.243 

170.243 

171.243 

8.83792 

5 89543 

0.82986 

0.75 

176.015 

174.015 

175.015 

7.47824 

6.6'3607 

0.77142 

O.BO 

179.600 

177.600 

178.600 

6.11856 

7.42270 

0.6^184 

0.B5 

182.854 

180.854 

181.854 

4.75888 

8.04565 

0.56406 

0.90 

185.644 

183.644 

184.644 

3.3992c 

8.53830 

0.42124 

0.95 

187.850 

185.850 

186.850 

2.0352 

3.87725 

0. 25666 

I .00 

189.368 

187.368 

188.368 

0.67984 

9.04:39 

0.07362 
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single  :c::tcn  bridge  considering  lass  of  lo^'d 


SINGLE  ;ONTCN  BRIDGE 


AirLI  D LOAD- Till;  CURVl.  FOR  H=Tq 


Equation  of  Motion: 

= (l+;3o<  (l+y3oL  ^)V^At 

.77  2^2  p 

‘ K' 


or 

Ayp 

=-■  By3_  + CVj_^t  + D 

Time 

1 

it 

) 

? = ^/20 

To 

A 

p 

C 

D 

7 

V^t 

0 

0 

1.025 

0.975 

1.000 

0 

0 

0 

0.05 

0.1 

1.055 

1.005 

1.030 

0.0017 

0.0016 

0.0032 

0.10 

0.2 

1.145 

1.095 

1.120 

0.0051 

0.0918 

0.0U8 

0.15 

0.3 

1.295 

1.245 

1.270 

0.0085 

0.0270 

0.0238 

0.20 

0.4 

1.505 

1.455 

1.480 

0.0119 

0.0574 

0.0370 

0.25 

0.5 

1.775 

1.725 

1.750 

0.0153 

0.1008 

0.0499 

0.30 

0.6 

2.105 

2.055 

2.080 

0.0187 

0.1566 

0.0617 

0.35 

0.7 

2.495 

2.445 

2.470 

0.0331 

0.0221 

0.0719 

0.40 

0.8 

2.945 

2.895 

2.920 

0.0255 

0.2997 

0.0805 

0.45 

0.9 

3.455 

3.405 

3.430 

0.0289 

0.3836 

0.0874 

0.50 

1.0 

4.025 

3.975 

4.000 

0.0323 

0.4738 

0.0929 

0.55 

0.9 

3.455 

3.405 

3.430 

0.0323 

0.5685 

0.0965 

0.60 

0.3 

2.945 

2.895 

2.920 

0.0289 

0.6643 

0.0952 

0.65 

0.7 

2.495 

2.445 

2.470 

0.0255 

0.7554 

0.0870 

0.70 

0.6 

2.105 

2.055 

2.080 

0.0221 

0.8340 

0.0700 

0.75 

0.5 

1.775 

1.725 

1.750 

0.0187 

0.3900 

0.0421 

0.80 

0.4 

1.505 

1.455 

1.480 

0.0153 

0.9121 

0.0019 

0.85 

0.3 

1.295 

1.245 

1.270 

0.0119 

0.8879 

-0.0500 

0.90 

0.2 

1.145 

1.^'95 

1.1.20 

0.0085 

0.8075 

-0.1107 

0.95 

0.1 

1.055 

] .005 

1.030 

0.0051 

0.6660 

-0 . 1723 

1.00 

0 

1.025 

0.975 

1.000 

0.0017 

0.4671 

-0.2256 
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appendix  B 


Determination  of  Damping  Coefficienta 
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Consider  the  single  iriass  oscillator  sho\"m  in  Fig.  3.  If  in 
addition  to  the  resistance  of  the  spring  the  iiass  encounters  a 
resistance  proportional  to  velocity  the  equation  of  motion  for  free 
vibration  (P=0)  is 

= Ky  + Gy  (77) 

The  solution  of  this  equation  is  of  the  form^^^ 

_Gt 

y = e ^ (A  cos  pt  + B sin  pt)  (78) 

where  A,  B and  p are  constants. 

The  solution  of  equation  (78),  for  a single  mass  oscillator 
subjected  to  an  initial  di^lacement  is  shown  schematically  in 
Fig.  6la.  The  equation  of  the  dotted  curve  connecting  the  maximum 
ordinate  is 

_Ct 

y = yo  ® ^ (79) 

or,  taking  the  natural  logarithm  of  each  side  of  equation  (79), 

C = -|!L  In  ^0/y  (80) 

Equation  (80)  can  be  used  for  the  experimental  determination  of  the 
coefficient  of  damping  C.  It  is  only  necessary  to  determine  by 
experiment  the  ratio  of  the  amplitudes  of  vibration  after  a period 
of  time  t. 

Data  from  the  tests  conducted  to  determine  the  virtual  mass 
coefficients  were  used  to  determine  the  damping  coefficients  for  the 
model  pontons.  An  example  of  such  data  is  shown  in  Fig.  6lb.  This 
trace  is  for  the  triangular  model  ponton  vibrating  at  a disilacement 
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of  2.06  in.  Considerinj  the  two  arbitrary  ;.oint.-  a ^nd  b,  0.2 
seconds  apart,  Uic  ratio  of  their  at:i[.litudes  is 


iz*o 


1.75 


The  tnacs  of  the  system,  including  the  mass  of  the  ponton  plus  one- 
half  the  ^ass  of  the  test  beam,  is  0.0479  Ibs-sec^/in.  whence. 


C 


to-ow) 

0.2 


or 

C = 0.26B  Ibs-sec^/in. 

Following  a similar  procedure  with  the  ponton  vibrating  in  air  a 
value  for  C of  0.0301bs-sec^/in.  is  obtained.  The  difference  between 
these  two  values  of  C is  the  damping  due  to  the  water  in  contact  with 
the  ponton. 

Similar  tests  were  conducted  on  the  rectangular  and  semi-circular 
pontons  at  different  displacements.  In  all  cases  the  damping  varied 
with  displacement  but,  to  simplify  calculations,  an  average  vtilue 
was  used  for  com.putations. 


I 
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APPENDIX  C 


Electronic  Analog  Computors 
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General 

An  electronic  analog  computor  deals  vd.th  problems  that  involve 
the  behavior  of  several  variables.  These  variables  are  represented 
in  the  computor  by  d-c  voltages,  which  are  referred  to  as  machine 
variables,  and  are  usually  related  to  the  original  variables  on 
some  convenient  scale.  The  given  mathematical  relations  between  the 
original  variables  are  expressed  by  an  analagous  set  of  relations 
between  the  machine  variables.  In  order  to  deal  with  a wide  range  of 
problems  the  mathematical  operations  on  the  machine  variables  are 
performed  by  basic  computing  elements  in  a small  number  of  simple 
operations.  In  general  a wide  variety  of  problems  can  be  treated 
withtthe  following  computing  elements: 

1.  Elements  which  will  multiply  a variable  by  a constant 
coefficient. 

2.  Elements  which  will  add  variables. 

3.  Elements  which  will  multiply  two  variables. 

4.  Elements  which  will  integrate  with  respect  to  time. 

Basic  Computing  Elements.-  The  addition  of  several  machine 

variables  V can  be  performed  with  the  network  shown  in  Fig.  62a.  Let 
A be  the  amplifier  forward  gain.  The  voltage  at  the  amplifier  input 
grid  is 

VVa 

From  Kirchoff's  first  law  the  sum  of  all  current  flowing  to  the  grid 
must  be  zero  or 


— 60" 


and  solving  for  Vq 


Vo 


(11  t 12  » 1>)  Aft. 

^1  ^2  ^3  (1.  + L +i.) 

® R2 


With  sufficiently  hiigh  amplifier  gain 


A 1 

and  equation  (81)  becomes 


Vo 


or 


where 


Vo  = -(^Vj^+  ^2* 

n = Ro/Rl  as  = R0/R2 


a-j  = HyR3 


(81) 


Thus  the  variables  are  not  only  summed  but  are  multiplied,  if  necessary, 
by  a coefficient. 

Multiplication  by  a constant  can  be  performed  with  the  network 
shown  in  Fig.  62b.  The  voltage  at  the  input  grid  of  the  amplifier  is 
Vq/A,  Equating  the  current  intc  the  grid  to  the  current  leaving  the 
grid 


(1°  - »o)l  * (?  - Vl)|-  - 0 

A Rq  A R|j_ 


and 


Vo  = Vi- 


(1  + A)^  + 1 


R. 


Since  A 1 we  have  to  a good  approximation 


Ro 

V = Vi  = a,V 


11 


(82) 


(83) 
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Fi^j.  62c  shov.'c  the  nctv;ork  for  multiplying  two  variables  tOjj;>ithei’, 
The  voltage  for  a setting  a of  the  voltage  divider,  is 

^o  = ^^1 

If  the  setting  a is  mde  proportional  to  some  voltage  by  means  of 
a servomechanism  then  the  output  voltage  aV^  vd.ll  be  proportional  to 
the  product  of  and  V2. 

Vo  = Vif(V2) 

Electrical  integrators  are  based  on  the  property  of  capacitors. 

The  current  through  a capacitor  is  proportional  to  the  derivative  of 
the  voltage  E 

i = CflE/dt  = CFE 

where 

P = d /dt 

Fig.  62d  shows  the  network  for  an  electrical  integrator  using  a capacitor 
and  a high  gain  amplifier.  The  equation  for  the  currents  at  the  input 
grid  is 

Vo  1 Vq 
(_  , (_  _V,)CP  = 0 

or 

V = ^ Vl 

° (i-A)ncp  + 1 

For  a large  aurplifier  gain  A 

~ RCP 

v/hich  corresi'onds  to  integration. 

Fig.  62e  shov/3  the  network  of  a sunriming  integrator  which  is  capaMo 


of  suirjiling  as  well  as  integration.  Folloivang  a similiir  derivation 
as  above  the  output  voltage  is 

1 Vp  Vo 

y . _ L (J-  + _f  + 

° CP  P-i  P-2  R3  (35) 

plock  Dia‘;rams.-  In  "setting  up"  a computor  to  solve  a mathematical 
equation  it  is  convenient  to  represent  the  relation  between  the  computor 
components  and  voltages  by  means  of  a block  diagram.  The  block  dia- 
gram describes  the  relations  between  voltages  and  computing  components 
in  the  sajiie  way  as  the  mathematical  equation  describes  relations  betv;een 
operators  and  variables.  The  symbols  commonly  used  to  represent  the 
basic  computing  elements  are  shovm  in  Fig.  62, 

As  an  example  consider  the  equation 

+ aY^  (86) 

Yj^  and  Y^  are  voltages  and  Yl  is  a function  ofY2.  The  voltage  Y2  v,’ill 
be  the  input  ''■oltage  to  a block  of  computing  components  as  shown  in 

Fig.  63. 


YiO 


<3  Y 


2 


FIG.  63  BLOCK  REPRESENTATION  OF  RELATIONSHIP  BET/JEEh' 
Y3_  AI4D  Y2. 

If  the  correct  combination  of  computing  elements  is  used  Y2^  will 


appear  as  the  output  voltage. 

Frequently  it  is  not  possible  to  solve  explicitly  for  a variable 
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in  equations  like 

Y^=  y2  + bY2  + cY^  (87) 

However,  since  appears  on  the  right  of  equation  (8?)  it  can  be  treat- 
ed as  an  input  voltage  to  the  block  of  computing  elements.  Since  Y-j^ 
at  the  input  terminal  must  equal  Y^  at  the  output  terminal  the  out- 
put voltage  is  treated  as  if  it  were  known  and  is  "fed  back"  to  the 
input  terminal  as  shown  in  Fig.  64. 


FIG.  64  BLOCK  REPRESENTATION  OF  RELATIONSHIP 
BETWEEN  Yj_,  Y2,  AND  Yo 

The  correct  relationship  will  be  established  if  suitable  computing 
components  are  used. 

Application  To  Floating  Bridges 

Equation  (5)  for  a single  ponton  bridge  considering  damping  is 

+ bdy^/dt^  + dy^  - f(t)  = 0 
or 

P^y^  + bPy^  + dy^  - f(t)  = 0 (88) 

where 

b = C/rfi^ 
d = (a22^+K)/mi 
f(t)  = F^(t)/m^ 

Solving  for  the  highest  derivative  of  y^  equation  (88)  bcco:;-os 
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P 7l  = - d7i  + - f(t)  (89) 

The  block  diagram  for  the  solution  of  equation  (89)  is  shovm  in  fig  65, 

The  voltages  bPyj^,  and  f(t)  are  applied  to  the  input  terminals 

2 

of  a summing  circuit  (Fig  62  a)  which  yields  - P yq.  These  input 

voltages^  however,  must  be  obtained  from  P^j_  by  integration  and 

2 

multiplication  by  suitable  constants.  The  voltage  - P y]^  is  con- 
sidered known  and  is  applied  to  the  input  of  an  integrater  which  yields 
Py]^.  This  voltage  is  multiplied  by  a constant  b (see  Fig  ^ b)  and 
serves  as  the  input  voltage  bPy^^  applied  to  the  summing  circuit.  The 

voltage  Pyj^  is  also  applied  to  the  input  of  another  integrator  which 
yields  - y^^.  The  voltage  - y^  is  applied  to  a phase  inverting  amplifier 
which  yields  y^.  This  voltage  is  multiplied  by  a constant  d to  pro- 
vide the  voltage  dy^^  which  must  be  applied  to  the  summing  circuit  and 
is  also  recorded.  This  record  constitutes  a solution  to  the  problem. 

The  voltage  -F(t)  applied  to  the  summing  circuit  is  supplied  by  a 
suitable  function  generator.  With  the  computor  set  up  in  as  shov.-n  in 
Fig  65  the  machine  variable  y^  is  forced  to  vary  in  the  manner  pres- 
cribed by  equation  (89). 

Ecjuation  (3)  for  a single  ponton  bridge  considering  virtual  mass 
can  be  written 


where  b = /fq 

8 = t k) 

f(t)  = F(t)  /M 


Solving  for  r equations  (90)  becomes 


p2y^  - dy^  + f(t) 


(91) 


Fig  66  shows  the  block  diagram  of  the  computer  set  up  for  the  solution 
of  equation  (91).  The  solution  is  carried  out  in  much  the  same  way 
as  that  of  equation  (88)  with  the  exception  of  the  addition  of  the 
servomechanism.  This  is  necessary  to  generate  the  product  of  by^^  and 

~ 8 y^  to  yield  the  voltage  — by^^P  which  is  applied  to  the  surniiing 


circuit. 

The  construction  of  block  diagrams  for  bridge  of  more  than  one 
ponton  presents  no  difficulties.  Equations  (2)  for  a two  ponton 
bridge  (n=2)  becomes 


^1  ^ ^11  yi  + ^2  ^2  - 8i(t)  = 0 
2 

^1  " ^21  yp  + b22  y2  - f2(t)  = 0 


(92) 


where  b-j^p  = 

(app  + k)  /Mp 

^12  " 

ai2  /Mp 

II 

X3 

^21  A‘i 

^22  " 

(^22  + ><)  Alp 

Fig  6?  shows  the 

computer  set  up 

for 

A separate  circuit 

is  constructed 

for 

to  the  summing  circuit  of  the  second  equation  is  taken  fron.  the  out- 
put of  the  first  equation  and  the  input  b^2  V2  first  equation  is 

obtained  fi'om  the  output  of  the  second  equation. 
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FIG.  5 nyOYAf'I  FORCE— DISPLACEMENT  RELATIONSHIP 

FOR  AN  M-4  PONTON 


SINGLE  PONTON  BRIDGE 
NEGLECTING  VIRTUAL  MASS, 
DAMPING  AND  MASS  OF  LOAD 


displacement— TIME  RELATIONSHIP  FCR  A SINGLE  PONTCN  BRIDGE 


FIG.  8 VIRTUAL  MASS— D I 


RATIO  OF  TIME  OF  CROSSING  TO  FUNDAMENTAL  PERIOD 

FIG.  10  RELATIONSHIP  BET.vEEN  IMPACT  FACTOR  AND  TIME  OF  CROSSING  FOR  SINGLE  PONTON 
BRIDGE 


RELATIGNSHIP  BETViEEN  IMPACT  FACTOR  AND  TIME  OF  CROScl^G  FOR  SINGLE  PONTCN  3RICGE 
CONS  I DERING  MASS  OF  THE  MOVING  LOAD 
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TIME,  MULTIPLES  OF  T^ 

DISPLACEMENT— TIME  RELATIONSHIP  FOR  A TV«  PONTON  BRIDGE 


TIME,  MULTIPLES  OF  T^ 

FIG.  14  DISPLACEMENT—TIME  RELATIONSHIP  FOR  THREE  PONT.N  3RIDG 


T-iREi  PONTON  aRIDGE  CONSIDERING 
VIRTUAL  MASS  AND  NEGLECTING 


T1^'E  RELATIONSHIP  FOR  THREE  PONTON  BRIO 
IRTUAL  MASS 


THREE  PONTON  BRIOC.E 
CONSIDERING  VIRTUAL  MASS 


TIME,  MULTIPLES  OF  Tq 

FIG.  16  DISPLACEMENT— TIME  RELATIONSHIP  FCR  THREE  PONTuN 

VIRTUAL  MASS  AND  DAMPING 
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TirE,  MULTIPLES  OF 

DISPLACEMENT— TIME  RELATIONSHIP  FCR  SINGLE  PONTON  BRIDGE  CONSIDERING  MaSS  OF  TnE  LOAD 


TIME,  MULTIPLES  OF  Tq 

DISPLACEMENT— TIME  RELATIONSHIP  FOR  SINGLE  PONTON  BRIDGE  CONSIDERING 
MASS  OF  THE  LOAD 
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OISPLACEMENT—TIME  RELATIONSHIP  FOR  T//0  PONTON  BRIDGE 
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TIME,  MULTIPLES  OF  T^ 

FIG  25  DISPLACEMENT— TIME  RELATIONSHIP  FOR  TJO  PONTON  BRIDGE  CONSIDERING 

VIRTUAL  MASS 


FIG.  25a  relationship  BEP^EEN  IMPACT  FACTOR  AND  TIME  CF  CROSSING  FCR  T'a'O 
PONTON  BRIDGE 
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FIG  26  DISPLACEMENT— T I Ml  RELATIONSHIP  t-OR  TWO  PONTON  BRIDGE  CONSIDERING 

virtual  mass  and  damping,  C = 100  LB-SEC/lN. 


TWO  PONTON  BRIDGE  CONSIDERING 
VIRTUAL  MAGS  AND  DAMPING  AND 
NEGLECTING  MASS  OF  THE  LOAD 
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TIME,  MULTIPLES  OF  Tq 

DISPLACEMENT— TIME  RELATIONSHIP  FOR  TV.O  PONTON  BRIDGE  CONSIDERING 
VIRTUAL  MASS  AND  DAMPING,  C=  1000  LB-SEC/ 1 N . 
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DISPLACEMENT— TIME  RELATIONSHIP  FOR  FIVE  PONTON  BRIDGE 
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COMPARISON  OF  NUMERIOAL,  EXACT,  AND  ANALOG  SOLUTIONS  FOR  A T\VO  PONTON  BRIDGE 


35  VIRTUAL  MASS  APPARATUS 
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37  TEST  PCraONS 
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Fia40  CALIBRATION  CURVE  OF  TEST  8EA1 


FIG.  47  RESULTS  OF  IMPULSE— D I SPLACEMEMT  TESTS  GU  TRIANGULAR  PONTOr 


RESULT  CF  IMPULSE— 01  SPLACEMEMT  TEST  ON  TRIANGULAR  PONTCN 
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RELATIONSHIP  BETvVEEN  IMPACT  FACTOR  AND  TIME  OF  CROSSING  FOR 
SINGLE  PONTCN  BRIDGE  CONSIDERING  MASS  OF  THE  LOAD 
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RATIO  OF  K TO  EI/L^ 

FIG. 57  RE'.ATIONSHIP  BETWEEN  INTACT  FACTOR  AND  RATIO  OF  K TO  EI/L^ 
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BLOCK  REPRESENTATIUR 


(a)  addition 


(b)  multiplication  by  a constant 


(c)  multiplication  of  two  variables 


FIG. 02  BAS  1 C COMPUT I NG  COMPONE NTS 
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NOTICE:  WHEN  GOVERNMENT  OR  OTHER  DIL\ WINGS,  SPECIFICATIONS  OR  OTHER  DATA 
ARE  USED  FOR  ANY  PURPOSE  OTHER  THAN  IN  CONNECTION  WITH  A DEFINTfELY  RELATED 
GOVERNMENT  PROCUREMENT  OPERATION,  THE  U.  S.  GOVERNMENT  THEREBY  INCURS 
NO  RESPONSIBILITY,  NOR  ANY  OBLIGATION  WHATSOEVER;  AND  THE  FACT  THAT  THE 
GOVERNMENT  MAY  HAVE  FORMULATED,  FURNISHED,  OR  IN  ANY  WAY  SUPPLIED  THE 
SAID  DRAWINGS,  SPECIFICATIONS,  OR  OTHER  DATA  IS  NOT  TO  BE  REGARDED  BY 
IMPLICATION  OR  OTHERWISE  AS  iN  ANY  MANNER  LICENSING  THE  HOLDER  OR  ANY  OTHER 
PERSON  OR  CORPORATION,  OR  CONVEYING  ANY  RIGHTS  OR  PERMISSION  TO  MANUFACTURE, 
USE  OR  SELL  ANY  PAI  ENTED  INVENTION  THAT  MAY  IN  AIU  WAY  BE  RELATED  THERETO. 
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